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IN MEMORIAM H. B. DORGELO 


The death on 6th March 1961 of Professor Dr. H. B. 
Dorgelo, rector magnificus of the Technological University 
of Eindhoven, was a severe loss to the Dutch scientific 
world in general and to this journal in particular. Dorgelo 
was one of the founders of Applied Scientific Research 
directly after the second world war, and as chairman of 
the Governing Board he actually helped to overcome the 
numerous obstacles that lie in the way of a new publication. 
From the very beginning the emphasis was laid on the 
international character of the journal and it may be said 
without hesitation that by now this intention has been 
fully realized. As General Editor I had often occasion to 
consult Dorgelo on whose full support I could always 
count. Now that I am taking over his task in the Governing 
Board and turning over the Editorship I wish to express 
the hope that Applied Scientific Research may be a lasting 
memorial to a man who gave himself wholeheartedly to 
his task. 


Rk. Kronig 
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HIGH-FREQUENCY DIFFRACTION BY A WEDGE 
WITH A LINEARLY VARYING SURFACE 
IMPEDANCE 


by L. B. FELSEN 


Office of Naval Research, London Branch, London, England 


Summary 


A previously derived formal solution for the two-dimensional Green’s 
function for a wedge with a linearly varying surface impedance is examined 
in the quasi-optic range of frequencies. An asymptotic representation, 
valid both in the illuminated and shadow regions, is derived for the scattered 
field due to an incident plane wave, and its properties are analyzed as a 
function of the rate of variation of surface impedance. A geometric-optical 
interpretation of the solution is emphasized. 


§ 1. Introduction. As noted in a previous publication 1), the 
two-dimensional problem of diffraction by a wedge with a linearly 
varying surface impedance (or admittance, depending on the 
polarization of the incident wave) can be solved by the conventional 
method of separation of variables, and the solution can be presented 
in alternative forms which are advantageous for different investiga- 
tions. A radial transmission formulation has been found useful 
for the analysis of radiation from a tapered surface wave antenna 2), 
the surface impedance of which varies in the above-described 
manner. In this paper, an angular transmission formulation is 
utilized to determine the scattering properties of the variable 
impedance wedge in the quasi-optic wavelength range. The dif- 
fracted field due to an incident plane wave, observed far from the 
edge, is expressed in the shadow region in terms of an asymptotic 
expansion involving terms of the form (kp)~”~!In(d, ¢’), m= 
=0,1,2,..., where k is the free-space wavenumber, (p, ¢) the 
cylindrical coordinate of the observation point and ¢’ the incident 
wave direction. The angular function I(¢, ¢’), which also depends 
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on the rate of impedance variation on the wedge, is given in terms 
of an integral and, for m > 1, is derivable via a recursive procedure 
from a knowledge of Jo(¢, ¢’). Explicit expressions for Ip are 
obtained for two special cases: 1) a wedge having a linearly varying 
surface impedance on one side and a “‘perfect absorber’ of the 
Sommerfeld type on the other side, and 2) a wedge having a purely 
reactive impedance with the same rate of (linear) variation on 
both sides, but with the algebraic sign of the reactance on one side 
being the opposite of that on the other. In the illuminated region, 
the asymptotic evaluation, carried out to O(1/,/kp), yields a 
reflected geometric-optical field in addition to the diffracted 
contribution. 

While the variable impedance wedge configuration described 
herein might conceivably serve as a model for the analysis of 
diffraction by sharp mountain ridges with certain variable proper- 
ties, or by high-speed vehicles in the upper atmosphere surrounded 
by a plasma sheath with a certain variable density, the main 
purpose of this investigation has been the comparison of the 
rigourous asymptotic solution in the quasi-optic range with pre- 
dictions based on generalized geometric-optical procedures. Such 
procedures have been employed with considerable success by 
Keller and his co-workers 3) for the analysis of diffraction by 
relatively arbitrary structures. In Keller’s geometrical theory of 
diffraction, the scattered field as k ~ co is viewed as arising from 
various types of reflected and diffracted rays, whose amplitudes 
are determined by the properties of the obstacle surface “‘at’’ the 
point of emergence of a given ray. If this model is to be applied to 
the analysis of diffraction problems at large but finite values of R, 
it is evident that the word ‘‘at’’ employed above should be replaced 
by “in the neighbourhood of’, with the expectation that a slow 
variation of surface properties in an interval of a wavelength does 
not invalidate the local analysis of the quasi-optic scattering 
process. Since there is no precise estimate of what constitutes 
a slow variation, it is pertinent to investigate specific structures 
exhibiting different types of variability of surface properties and 
to examine the dependence of these solutions on the rate of variation. 
The present problem is suited to an analysis of the effect of a variable 
surface impedance in the vicinity of an edge singularity. 

If the surface admittance varies like Ap, where A is a constant 
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(which may have different values on the two wedge faces) and p 
is the distance from the apex, it is found that the diffracted wave 
amplitude for A <1 is essentially the same as that for a zero 
admittance wedge (although the wedge admittance tends to 
infinity as p > co), while the diffracted wave amplitude for A > 1 
is essentially that for an infinite admittance wedge. Ad <1 and 
A®*> 1 correspond, respectively, to slow and rapid admittance 
variations near p= 0. One notes that in the former case the 
diffracted wave amplitude is determined properly by the zero 
admittance value “‘at’’ the point of diffraction p = 0, while in the 
latter case the edge admittance is of no consequence. For the 
two special cases mentioned at the end of the first paragraph, 
a study of the diffracted wave amplitude as a function of A can 
be carried out explicitly to determine over what range the 
quasi-optic approximation (involving the point of diffraction 
only) yields results acceptable to within some definite criterion 
of accuracy. 

The utility of an angular transmission analysis of certain wedge 
problems, involving the use of the Kontorovitch-Lebedev 
transform, has been emphasized previously 4) ®). In the present 
context, the angular domain is not restricted to the region 
0 <q <« bounded by the two half planes forming the wedge, 
but is allowed to range from ¢ = — co to ¢ = oo. The requisite 
boundary conditions at ¢ = 0, « are then restored by superposing 
in the infinite angular domain an infinite set of ‘‘image’’ solutions. 
While the infinite angular space concept has been exploited 
primarily for the analysis of diffraction by curved objects §) 7), 
its considerable utility in connection with structures possessing 
edge or tip singularities 1) is not to be overlooked. In fact, a prime 
virtue of this type of analysis (not sufficiently emphasized in 
the literature) is a resulting formulation which permits the 
immediate identification of geometrically reflected, and of dif- 
fracted, field contributions in an asymptotic evaluation at high 
frequencies. 


§ 2. Formal solution. It was shown in 1) that the electro- 
magnetic fields radiated by a line source of electric or magnetic 
currents in the presence of the variable impedance wedge (see fig. 1) 
can be inferred from a scalar Green’s function G. The impedances 
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0,0 at the wedge faces d = 0, « have the form 


ea Zo tkpco,x for a line source of electric currents, 
OS ia 5) : . 
Zo/tkpco,x for a line source of magnetic currents, 


where Zo = 1/uje is the characteristic impedance of free space, 
p the radial variable, k the free-space wavenumber, and co, cy are 
constants restricted such that Imco,, <0. (A time dependence 
exp (— zwt) is assumed throughout). G satisfies the boundary 
condition G = + ¢o,4(@G/éd) at 6 = 0, a. 


REFLECTED RAY ae 
BOUNDARY 


p-a 
Fig. 1. Physical configuration and various geometric-optical domains. 


For a study of the high-frequency diffraction properties of the 
wedge surface, the Green’s function G is most conveniently expressed 
in terms of the angular transmission (radial eigenfunction) re- 
presentation in (22) of 1): 

G(p, p’) = 


ico 


=f fp(1—e™") Hy) (hp) Hal? (Rp') Galo, $512) du, |b — $'| > 2. (1) 
0 


The angular characteristic Green’s function Gy is given in closed 
form by (12) of 1): 


(e##< os eae tneres) me t's eiMla—o>)) 


ea = Qine™ — Tone) - 
Co — | tly — | 

= = —_— _, Ixy = l(u) = ———__> 3 

Io = To(u) eee (1) Tyee (3) 


where ¢— and ¢, denote the lesser and greater, respectively, of 
the variables ¢ and ¢’. Alternatively, Gg can be expressed in terms 
of images in an infinitely extended angular transmission line as 
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imei )poteuls 
Cpl, Bs H?) = Cpoolh ) + ue 0” TDa"G goo(h, 2ne — P') > 


co 


+ 2 (Lalo) "1G geolP, 2ne + $') + Cool, —2ne +O] 


v=1 


+> Peele Gee: —2na — ¢'), Inmu> 0, (4) 
n=0 
where 
eitld—o'| 
Grol, $') = —, Imp > 0. (5) 


The physical interpretation of (4) has been discussed previously 
(see fig. 4 of 1)). The representation in (1) is restricted to the 
angular region |é — ¢'| >a, 1.e., to the geometrical shadow 
region of fig. 1 from which the source is not directly visible. 


$3. Asymptotic evaluation. 

1. In the shadow region. When p’ + oo (plane wave inci- 
dence) and kp > 1, we may represent the Hankel functions in (1) 
by their large argument asymptotic formulae 


Hy (kp’) ~]/ 


2 


eilke’ —tyun— jn) 
mtkp’ 


p' > ©, (6) 


titers es) (u, m) 
H,Y(k ~) i(ke— tun — tn) ek ee es 
u (Rp) ae = (— 2ikpym p> |ul, (2) 


_ 4p? — 19)[4u? — 37]...[4a2 — (2m —/1)?] 
_ 22mm | 
Since the integrand in (1) is exponentially damped for large values 
of « when |¢ — ¢’| > a (see (9) of 1) and (4)), the error incurred 
by employing (7) over the entire range in 4 can be shown to be 
exponentially small and therefore negligible. The field in the 
shadow region is therefore given asymptotically by 

Inl¢, $') 


(é : t muiliA kh , = eke OR ere O yews 
(p, p’) ~ A( eee, oe en kp> 1, (9) 


(tu, m) (u, 0) =1. (8) 


where 


aoe 
Ay = V2 ete, (10 
TCX 


In(d 8) = — 46 f mle m) sin ual, 8s 2) Qu. (11) 
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Since 
l 
(i, m) = “5 [u? — 4 — m(m — 1)] (u, m — 1), (12) 
and (d?/dg*)Gy = — wG,, ¢ #¢’, it follows that the I, can also 
be obtained from Jo through use of the recursion formula 8) 
l 2 
In = m oe + m(m — 1) + | dehy et = | ee aS) 


From the nature of the factors A(kp’)A (Rp) in (9) it is noted that 
the diffracted field in the shadow region has the form of an outgoing 
cylindrical wave, and can be viewed as arising from the incident 
ray which strikes the edge and is scattered radially in all directions 
as shown in fig. 2. The amplitude of the diffracted rays as a function 


as Rays 
kon /o\ 
Eee ala 


Fig. 2. Quasi-optic field contributions. 


GEOMETRICALLY 
REFLFECTED RAY 


of ¢, ¢’ and the surface properties is given by Im/(¢, 4’). Of particular 
interest is the dependence on the surface impedance 4,4. Although 
the diffracted rays appear to emanate from the edge, their amplitude 
is determined by the impedance in a neighbourhood surrounding 
the edge. This becomes evident from an examination of J, for the 
two extreme cases ¢o,, > 1 and co,x <1. For the case of a magnetic 
line source, for example, 49,4 co as p > O and 29,1 + 0asp — 00, 
for any finite value of co,4. When co, < 1, the impedance decreases 
very slowly away from the edge; in the evaluation of Im, one may 
set 0,4 * — 1 since the major contribution from the exponentially 
decreasing integrand arises from the range of small w. The resulting 
diffracted wave amplitude is approximately that appropriate to 
an infinite impedance wedge. However, when co,, > 1, the impe- 
dance decreases very rapidly away from the edge and the resulting 
Im is approximately that corresponding to I0,q ~ 1, i.e., a zero 
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impedance wedge. Thus, it is noted that for slow impedance 
variation (per wavelength) near the edge the diffracted wave 
amplitude is essentially determined by the impedance “‘at’’ the 
edge, while for a rapidly varying impedance the properties at the 
edge are of no consequence and the diffracted wave amplitude 
is determined by the properties of the surface away from the 
edge. 

To obtain a more precise estimate of what constitutes a slow or 
rapid impedance variation as employed above, we consider an 
auxiliaiy problem for which I’, = 0. This implies that the surface 
at ¢ = a, instead of having a surface impedance Z,, is a “perfect 
absorber for angularly propagating waves’’ (see remarks preceding 
(21) of 1)). Gg as represented in (4) now comprises only the contri- 
butions from the source at ¢ = ¢’ and an image source at ¢ = — @’. 
Upon writing 


2 
| 14 
P aco + 1 a 


one can evaluate Jo in (11) as follows: 


To(p, $') = 


2, ae Si bad pedis ] 
alz+b—#1 a1 +18) 2—G+4) 
1 
—4-|BG+$' 200) BG +4 +20) ], Pa=0, (15) 


where 


rors 50 


a e ” 
= ete F (is 2) Ex(y) = 
i(——-), Ex(y) i} —dy. (156) 


The exponential integral Ej(y) is a regular function of y in the 
complex y-plane cut along the negative real y axis. (150) can likewise 
be considered as the analytic continuation of (15a) for complex 
and negative real values of x. E1(y) possesses the series represen- 
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CO 1) Fate 
Ex(y) el y In NY Ss ( y) 


n-1 nin 


-, y=0.5772. (16a) 


For large y, successive integration by parts in (15d) yields the 
asymptotic formula 


og oes m! : 3m 
Cae erga = (— ym” Pea ae er aL Oe 


The exponential integral can also be expressed in terms of the 
sine and cosine integrals 


y y 


sin ¢ St 
Si (y) = | —_ Pee Ge Hy ts | — dt, (16c) 
0 fore) 
as 
Ex(— iy) = —Ci(y) — 1 Si y) + x (164) 


For large values of co, use of (16a) yields 


Ban SSeS — ic E a L-() ( ~ ) ,co> 1, (17a) 


Aco XI 


while for small values of co (160) leads to 


] 
= 3: [B(x1, co) — B(x2, co)] = ACs _) + O(co), co<1. (176) 

A plot of Q(f) = B e-8E;(— 8) for positive real, negative real, 
and negative imaginary values of 6 is shown in fig. 3, where 
B = (x/co). Since x stands for (6 + ¢’ + 2) and is always real, 
the range of real values of # allows for positive and negative real 
values of co, i.e., purely reactive surface impedances. The negative 
imaginary values of f, on the other hand, account for positive 
real surface impedances when x > 0. It is noted for all cases that 
\O| < 0.05 for |B| < 0.01 and (1 — |Q|) < 0.05 for |f| > 20. To 
facilitate an interpretation, let us rewrite (15) as 


Brees a 1 | ) 
= Bae] za 


— 5 [1 + 20(6)] + Fag ll + 20lpa)), (184) 


(ps) =p « # Ei(—p), fis = <a ,m2=—¢4+4' Fa. (180) 
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It is difficult to predict without detailed calculation via (18a) the 
effect of various values of co on the magnitude of Jo since the 
pertinent parameter is 6 which contains both 4, ¢’ and co. How- 
ever, we can make the following statements regarding maximum 
effects:a) If |61,2|<.01, the coefficients of (8x1,2)—1 differ in magnitude 
by at most 10 percent from unity. A coefficient equal to unity 
represents a value co = oo, 1.e., a surface with infinite and zero 
impedance, respectively, for the electric and magnetic line source 
excitations. b) If |f1,2| >20, 20(61,2) differs in magnitude from the 
value (— 2) by at most 10 percent. A value of (— 1) for the co- 
efficients of (8x1,2)~! represents a value co =0, corresponding to a 
surface with zero and infinite impedance, respectively, for the 
electric and magnetic line source excitations. 


P=an 
Anns 
En 


(eS ee 


A 
A 
sa el 


Fig. 3. Plot of Q(f) = B e-8E1(— 8). 


2. In the illuminated region. To carry out an asymptotic 
field evaluation in the illuminated region |¢ — ¢’| <q, it is con- 
venient to separate from the image representation (4) those terms 
for which | — | <a, where @ is the location of a given source 
in the infinitely extended angular region. These terms alone can 
give rise to a geometric-optical field contribution, as shown below. - 
If the wedge angle « > z, the only sources in this category are 
those situated at p = + ¢’, 2« — ¢’, where 0 < (4, ¢’) < «. Thus, 
we write 


Ge = @, + (Gs — G), (19a) 


where 


Gy = Cools $') + PoGpoolh, ~$') + DaG goal, 2c — $'). (198) 
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Since (Gg — Gg) behaves for all 0 < (¢, $’) <a like exp[— (Imy)}, 
B > 2, when «too, an asymptotic evaluation of the contri- 
bution from this term to the integral (1) is carried out as in (9) 
and yields only a diffracted wave effect. For an evaluation of the 
contribution from Gs in (196), we examine separately each of the 
terms. 


Contribution from G,,,(¢, ¢’). Upon substituting G,,, into (1), 
replacing H,")(kp’) by its asymptotic form (6), Hy (kp) by its 
Sommerfeld integral representation 19) 

6 —ico 
] 
Hy (kp) = — fav Cc ery, joe UT, (20) 
A 
—6+ oo 
performing an interchange of the orders of integration (permissible 
if |p — ¢’| > a+ 6) and carrying out the integration over yp, 
one obtains for the integral L; defined as follows: 


Lal = $1) = EJ Gals ul — 0) Hyp) Hyp") du = 


= A(bp') = | dpe *D(y, | 4"), p! 00, |6—$'|> 4 +0. (21) 


—6+100 


1 

p= . (21a 

ee ete a Sore aye 
Unlike the left-hand side of (21), the integral on the right-hand 
side is defined as well when the restriction |¢ — ¢’| > x + 6 is 
removed. If the path of integration in the complex y-plane is 
chosen symmetrically with respect to the origin, the value of Ly 
for | — ¢'| <2 differs from that for |p — ¢’| > a by the residue 
at the pole yp = 2 — |¢ — ¢’| as shown in fig. 4. Thus, for any 


6-4! 
Lali — #1) =~ Alle’) | e-Hreone-#) Hl — 6 — $1) — 


6—ico 


u: fay eike cos w D(y, ld a | } (22) 


—6+%00 
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where 

FL (deeb ep cen pet ae ease (224) 
It is noted from (22) that G,,, contributes the incident field 
in the illuminated region |6 — ¢’| < a, and that the diffracted 
field as represented by the integral has the same functional form 
(as a function of |é — ¢’|) for any |f — ¢’|. (For an incident plane 
wave of unit amplitude, the normalization }7A (kp’) = 1 is employed). 


Yep 


@) lp-f'|> 1 (b) IP-P'l <7 


Fig. 4. Paths of integration in complex J-plane. 


A steepest descent evaluation (through the saddle point at y = 0) 
of the integral in (22) for large values of kp yields to O(1/+/Rp) 
and for all values of |6 — ¢’| 14) 


Lali — $')) ~ A lk) | em Hre-*) Hla — jp — $l) — 


” | ] 
MAS 
2 Cam sre Nheete 


)-re.i§— #0}, @) 


where 
Re ete ie | Fe) ein | ooh alan Se 
p> > 88 Red easy oll bsg ara) |, (23a) 
sgn (t# —y)=+1, Ex. (230) 
F(é) is the function 12) 
7 2 ize? 2? 
(&) = fe € dxe iy (23c) 
: (to 
which behaves for large & like 
tin 
PE) oe 


The first term inside the square brackets in (23) represents the 
incident field, the second term a diffracted contribution and the 


HIGH-FREQUENCY DIFFRACTION BY A WEDGE 181 


third term the transition effect near the shadow boundary 
lp — ¢’| =a. The extent of the transition region is defined 
approximately by € < 4; for > 4, the asymptotic formula (23d) 
applies with good accuracy so that T in (23a) vanishes in that 
range. Since the transition region at |¢ — ¢'| ~ a arises entirely 
from the contribution due to G,.,(¢, ¢’) and does not contain a 
dependence on the impedance properties of the wedge, one notes 
that the light-shadow transition behaviour is independent of the 
wedge impedance. 


Contribution from 1 G,,.(¢, —¢'). Upon employing (14) 
and following the same procedure as above one obtains for the 
integral L» defined in the range (¢ + ¢’) > a as 


L(¢+¢') = EJ ToC salts —$')u(1—e") yD (kp) Hy (kp’) du (24) 


the result (p’ — oo) 
Lop + $') =Lig+¢)+Le(d+¢); G+¢) > (25a) 


6 —ico 


d ; etke cos o 
La'(b + $') = — Ale’) -— | dpe 


—6+i00 


[Bd + ¢’ + p— 2, co) —B(6 + 6’ +p +2, 60)], (258) 


where Li(¢ + ¢’) and B(x, co) are defined in (21) and (15a), re- 
Pee The integrand in (255) possesses branch points at 

= +2 ~ (¢ + ¢’) in view of (15d) and (16a). When (6 + ~)> a2 
ot branch points occur to the left of the symmetrically chosen 
path of integration in the complex y-plane, and L’ is defined 


Vere CH +P? 


Fig. 5. Contours for J’ when (¢ + ¢/) <2 


as in (250). When (¢ + ¢’) <2, the value of L»' differs from 
that in (25) by the branch point, and associated branch cut, 
contribution due to yp =a — (+ ¢’) as shown in fig. 5. The 
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branch cut has been drawn so as to assure convergence of the 
integral along cp as determined by the behavior of exp(7kp cos ). 
Thus, for any value of (6 + ¢’), 


Lo(é + ¢’) = 
—11(¢+¢')+Le'(6+4’)+Le"(6+¢')A(a—(¢+¢4')), p’ +00, (26a) 
Lo"(¢ + $’) = 


l at ' 
= — Alte) — | ay oH * Bey — yo, co), yo=a— (6 +4), (260) 
Astco 
Co 
where L; and L’ are given in (22) and (250), respectively. 
To effect a simplification of L2” we introduce the change of 
variable 


y= Yo — 2 (27) 
so that 
_ 1A (kp’ 
Le"(¢ +d) = et fe: eike contte= 17) Bi — sz tca)t (28) 


Co’ 
where cy’ is a contour enclosing the branch cut along the positive 
real z axis in the positive sense. Upon substituting for B(— 72, co) 
the exponential integral from (15d) and then the series expansion 
(16a), one notes that the only contribution to c,’ arises from the 
(In z) term since all the other terms in the series are continuous 
across the positive real z axis. Thus, 


1tA(R ; Ses 
Tah teddy a fa: efiv costo) gizlee nz, (29) 
Finally, since ' 
In (x e”?") = In x + 12m, (30) 
one obtains for (29) 
2 A(kp’ ‘ : . 
Lab + $') = — 5°? [arexplitp coslye — ix) + ixleo] = (ta) 
0 


__ Alte! | x 
=— dx exp | 2 | kp cos pp cosh x + —— cosy} — 
2co ; |co| 
0 


: : % het 
- (ip sin pp sinh « + rT sin ») | , (316) 
0 
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where 
co= |cole"”, O<y <a. (31c) 


Since sin yp > O in the domain of existence of L2”, the integrand 
is exponentially damped with increasing ~. 

An asymptotic evaluation of Le in (26a) for large values of kp 
can now be carried out directly. The result for Li(f + ¢’) is given 
in (23), and a saddle point evaluation of L2’ in (25d) to O(1/4/p) 
yields (B is a slowly varying function compared with exp(ikp cos y)) 


Lo'(6+$") ~ 
1 
Joe Alte") (kp) (BYG+$! — 2, c0) — BUS + 4" +, c0)]. (32) 


_For an asymptotic evaluation of L2” one notes that the maximum 
contribution occurs in the neighborhood of x = 0. Upon expanding 
the exponent to O(«) and carrying out the integration one obtains 
(note: yy =a — (6 + 4) 


L2"(¢6 + ’) ~ — $A (Rp’) eH? C89) AMEE YG Oa (33) 
Thus, 
Lal-+48)~1 Ae) {Folhosin(d +g lemon“ lo —$)— 
l 2 


1 1 
Feed aaa gh oie suring al 


(Bb + 4! — 1,00) — BG +4 +200) |—T d + (4 


where J(u) is defined in (3). 

It is noted that for |f — ¢’| > a, (34) together with (23) agrees, 
as it must, with (15) and (9). Thus, the functional form of the 
diffracted wave contribution is the same for all geometric-optical 
domains. The expression for Lz in (34) becomes infinite as 
(¢ + d) +2, since the asymptotic evaluation of Le’ in (32) does 
not apply when the branch point approaches the saddle point at 
y = 0. However, the divergent term increases only logarithmically 
so that the affected wave amplitude behaves like 


(1//Rp) In (x — $ — $’). 
Since kp 1, the excluded range of (x — 4 —¢’), in which 
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| In (« —¢—¢')| is O(+/kp), is therefore very small. As in (9) 
and (23), the term associated with A(kp) in (34) represents the 
diffracted wave which appears to arise from the edge ; its dependence 
upon the wedge impedance function co has already been discussed 
in connection with (15). The transition function T(p, ¢ + ¢’) 
compensates for the divergent behavior of that portion of the 
diffracted wave amplitude which behaves like 1/(x — ¢ — ¢’) as 


(p+ $) +2. 


F sin(p+#') = CONST- 


INCIDENT 
REFLECTED RAY ag 


a aS 


p=a 


Fig. 6. Interpretation of geometrically reflected rays. 


The first term on the right-hand side of (34) exists only when 
& < (a — ¢’); the exponential term exp[— (ikp cos(d + ¢’))] per- 
mits its contribution to be identified as that from a ray which is 
reflected from the wedge face at d = O at the angle of incidence 
(see fig. 6). The amplitude of a reflected ray, given by 


icokp sin(? + ¢’) — 1 
icokp sin(d + ¢’) + 1 ’ 


is constant along the ray since the equation of a reflected ray is | 
given by psin(¢? + ¢’) = constant, as noted from fig. 6. The 
amplitude differs, however, for each ray so that the complex of 
reflected rays does not define a uniform plane wave. It is of interest. 
to compare Jo in (35) with the plane wave reflection coefficient | 
appropriate to an infinitely extended plane surface with constant 
impedance Zo, : 


Tolkp sin($ + $")] = (35) 


7 y sin ¢’ — | 
i A= —— == 
o(#) ysing’ + 1 


, (36) : 
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where the constant 7 is defined for the two polarizations as follows: 


ger: 
P Aieoeeg = [see 29, for E parallel to plane, (36a) 
Le 


l 
y = —., for H parallel to plane. 
20 


¢’ is the angle of incidence measured away from the surface. 
Along the surface ¢ = 0, (35) can be written as 

y(x) sind’ — 1 
y(x) sind’ + 1 
x is the distance measured from the apex to the point at which the 
specularly reflected ray passing through P(p, 4) leaves the surface. 
Upon comparing (36), (36a), and (37), one notes that the amplitude 
of the reflected ray is the same as would be obtained for a uniform 
plane wave reflected from an infinite surface having a constant 
impedance identical to that at x. This conclusion is in agreement 
with predictions based on geometric-optical reasoning. 


Io(kp sin ¢’) = a Vy () = theox, x = pcos. dé. (37) 


Contribution from 1G,,.(¢, 2x — ¢’). The contribution from 
the image source at (2% — ¢’) in the infinitely extended angular 
region is given by Lo(2« — ¢’ — ¢), where Lg is taken from (34), 
provided that co is replaced everywhere by cy. The physical inter- 
pretation of the various contributions is the same as above, except 
that all considerations now involve the wedge face at ¢ = a. 


§ 4. Final result. The asymptotic solution for all values of ¢ 
and ¢’ (except in the transition regions ¢ + ¢’ ~ a and 20 — 
— (d+ ¢’) ~ a) can now be written to O(1//kp) for p’ > co, 
ke => 1,/a'>" as: 

a / , / 
G(p,p') ~ 7 AlFp'Gope(p, $) +Gaitele, #') + Ctrans(e, #)1, (38) 


Gopt = 7 0°) H(z — |p — $'|) + 
+ Dolkp sin($ + $')] eM ees A (a, — $ — $') + 

+ Talkp sin(2«a—¢ — $’)] ei e8@2—9-# Hix, — (2a — $')], (38a) 

Gaise = — 4ilo(¢, $') A (Ap), (380) 

Girone = — Te, 6 — $'!) — Tle. 6 +’) — T(p, 2 —$ — 4’), (380) 
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where A(x) is defined in (10). The diffracted wave function 
To(¢, ¢’) has been shown to have the same functional form (as a 
function of ¢, ¢’) for all geometric-optical domains. An integral 
representation for Jg9 valid when |¢ — ¢’| > a is given in (11). 
Thus, Zo could be evaluated from (11) and the resulting function 
employed for all values of ¢, $’. Gipang NM (38c) is significant only 
in the various transition regions. It affords a complete transition 
through the light-shadow boundary |¢ — ¢’| =a, but does not 
apply in the immediate vicinity of the reflected wave boundaries 
at (6 + ¢’) =a, (24 — d — ¢’) = a (see remarks following (34)). 

One verifies readily that the asymptotic solution in (38) satisfies 
to O(1/4/kp) the required boundary condition G = F c¢o,.(éG/é¢) 
ato =O ose Since Gs satisfies this condition, one notes that for 
Io in (11) 


Io =F 2 
Ze car isos 
: Oe d¢ 
For the geometric-optical term in (38a) one finds when ¢’ < z, 


thule mn Got (1 co cot ¢’ 
Con ie dd icokp sin ¢’ + 1 


bs: (39) 


yo, $=0, 0 


and similarly for 6 = «. Since dG,,,/d¢ is 0(1) in kp, the geometric- 
optical contribution satisfies the boundary conditions at the wedge 
faces to O(1/kp). 

It is of interest to note that the asymptotic solution in (38) does 
not exhibit explicitly the presence of the surface wave discussed 
in 1), which can propagate along the wedge face at ¢ = O when 
co > 0, or along the wedge fact at 6 = « when cy > O. It is recalled, 
however (see (20) of 1)), that the surface wave varies like 1//kp 
as p -> co so that its effect in this case is indistinguishable from 
the diffracted field contribution. 

An explicit evaluation of Jg has been obtained in (15) for the 
special case [y= 0 where the wedge face at = « is a perfect 
absorber for “angularly propagating waves’. In this connection, 
we point out that Jo in (11) can also be evaluated for the physically 
more significant special case wherein the wedge faces at o=0 
and @ = « are purely reactive, with the reactance at 6 = 0 opposite 
in sign to that at $ = «, i.e., ¢ = —c, > 0. This implies via (3) 
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that Il’, = 1 so that from (4) 
ColPs Ps w?) = Gg(h, Os w)| rer. + 


] co co 

+ @ ~ ) Sy Gracsldaenar—d obs (T'y = 1) Sh Galt, —2na 4, 
0 n=1 n=0 

co> 0. (41) 


Thus, Gg can be expressed in this case in terms of the result for a 
perfectly reflecting wedge (Ig = Ix = 1) plus correction terms 
which involve only the first power of Ip or (1/Io). Upon sub- 
stituting (41) into (11) one obtains: 


if co 
Lo(d, ¢') = Lo(¢, $’)|n=r,-1 + eg. 2 Blane + 2+ 9, Co) — 


— B(2nz —2 + 9, co) + B[2(n + lla —x— , — co] — 
— Bi2(n + lla +x—, —co}}, p= 44+ 4, (42) 


where the first term on the right-hand side of (42) is the known 
expression for the perfectly reflecting wedge 13) 


Ae 2 
To(¢, P)\ro=re=1 eS Ba sin - : 
l * 
2 4 2 : ) , (42a) 
4 Jt X a 
cos —- — cos—(¢— ¢’) cos— —cos—(¢ + ¢’) 
a a a is 


_ and the remaining terms are corrections involving the exponential 
integral in (15a). 
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*) Substitution of (4), with 9 = Ix = 1, into (11) yields for Io a series representation 
which can be summed into the closed form (42a) via the formulas 
x-y 2 sin *% 


44 
+ cot 


1 co 
Coty = —— a 7 and cot 


% n=1 #2 — (nz)? cos v — cos x" 
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SOME DIFFRACTION PATTERNS OF AN 
ABSORBING HALF-PLANE *) 


by CHESTER J. MARCINKOWSKI 


Microwave Research Institute, Polytechnic Institute of Brooklyn, New York, U.S.A- 


Summary 


Far-field diffraction patterns are presented for the two-dimensional 
electromagnetic problem of plane wave scattering by a half-plane with a 
constant surface impedance. Previous results have been mainly concerned 
with a small surface impedance . In contrast with this, our diffraction 
patterns haye been computed for the relatively arbitrary impedance of a 
half-plane which acts like a black screen for a prescribed angle of plane wave 
incidence. Such a half-plane may produce no reflected plane wave. Numerical 
calculations for the far fields show that the back-scattered fields are more 
noticeably influenced by the surface impedance than the forward-scattered 
fields. Nulls are found in the back-scattered fields. This type of black screen 
produces approximately the same scattering as the Sommerfeld black screen. 
Simple qualitative estimates of the diffracted far fields may be made without 
elaborate calculations. 


§ 1. Introduction. The two-dimensional electromagnetic problem 
of diffraction by a dissipative half-plane with a constant surface 
impedance has been the subject of extensive investigations. 
Complete solutions for plane wave scattering were obtained by 
Senior 4) in 1952 by means of the Wiener-Hopf technique. These 
results were later extended by Senior 2) to solutions of the three- 
dimensional electromagnetic problem. It is well known that similar 
mathematical problems also occur in the coastal refraction of 
electromagnetic waves. A large amount of work has been carried 
out in this field particularly by the Russian school of investigators 
who include V. A. Fock) and his associates. They have worked out 
many details of this basic problem, including the Wiener-Hopf 

*) eranikcdacis partial fulfillment of the requirements for the Ph. D. degree at New 


York University. Sponsored by Air Force Cambridge Research Laboratories, ARDC, 
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solutions. The two-dimensional electromagnetic problem of. plane 
wave scattering by a constant impedance wedge of arbitrary angle 
was solved by Senior4) and Williams 5). For the special case 
of a wedge of zero angle the solutions of Senior and Williams 
simplify to those of the constant impedance half-plane. 

Previous investigators have appeared to consider their constant 
impedance, half-plane solutions primarily as approximations to 
a metallic half-plane with a small surface impedance. In contrast 
with this point of view we consider an absorbing half-plane with 
a relatively arbitrary, constant surface impedance. In particular 
we investigate some of the scattering properties of a half-plane 
which acts like one type of black screen with a constant surface 
impedance °). For a prescribed angle of plane wave incidence this 
type of half-plane produces no reflected plane wave. The three- 
dimensional electromagnetic problem of plane wave scattering by 
a black screen of this kind located on one side of a perfectly con- 
ducting right-angled wedge was solved by Karal and Karp’). 
Another type of black screen is a half-plane which is perfectly 
absorbing for “angularly propagating waves’’ incident upon it; 
this is Sommerfeld’s black screen 8). The surface impedance of 
such a half-plane is a function of position and its physical realizabili- 
ty is open to question. Nevertheless, the diffraction pattern of such 
a black screen is very simple to obtain 9) so that it is easy to 
compare the nature of the scattering produced by the two different 
kinds of black screens. 

It is characteristic of the more difficult Wiener-Hopf problems 
that they lead to factorized functions whose explicit evaluation 
requires the integration of a known function. For the constant 
impedance half-plane Senior 1) obtained an explicit evaluation of 
the integral in question. The Russian investigators studying problems 
in coastal refraction have made use of a much simpler solution 
obtained by Grinberg and Fock 8) with applications primarily 
to low impedance surfaces. We make use of the solution obtained 
by Grinberg and Fock and we extend the range of all parameters 
in their solution into the desired regions by analytic continuation. 
As a result of this procedure we obtain reasonably simple, explicit 
expressions for the two-dimensional electromagnetic problem of 
plane wave diffraction by a half-plane with a relatively arbitrary, 
constant impedance. In these expressions the surface impedance 


SOME DIFFRACTION PATTERNS OF AN ABSORBING HALF-PLANE 191 


is described by means of an “‘angle of perfect absorption”. This is 
defined as the angle at which an incident plane wave produces no 
reflected plane wave. This procedure has not only a simple geometri- 
cal significance but it also describes the half-plane impedance in a 
form which is useful and convenient for a physical interpretation 
of the diffraction patterns. Consequently, it is possible to describe 
the far-field diffraction patterns in terms of three simple dimension- 
less parameters: the angles of incidence, observation and perfect 
absorption. 

The diffracted far-field patterns for several representative angles 
of incidence and perfect absorption will be described. It was found 
that the back-scattered fields are appreciably influenced by the 
surface impedance; the forward-scattered fields are quite insensitive 
to the nature of the surface impedance, especially in the shadow 
region. The Sommerfeld black screen and the half-plane which 
produces no reflected plane wave (the “perfectly absorbing”’ half- 
plane) both give approximately the same diffraction patterns. An 
interesting feature is the existence of nulls which occur only in the 
back-scattered fields. The locus of these nulls is given by a simple 
equation. The results obtained allow us to interpret the asymmetric 
diffraction patterns for the absorbing half-plane as a kind of 
distortion of the symmetric diffraction patterns for a perfectly 
reflecting half-plane. In this distortion the locus of the nulls is 
shifted and the scattering amplitude is correspondingly distorted. 
The main properties of the diffracted far-fields may be predicted 
in a qualitative manner by means of relatively simple calculations. 


§ 2. The far-field solutions. In this section the asymptotic far-field 
solutions will be given for the two-dimensional electromagnetic 
problem of plane wave scattering by an absorbing half-plane with 
an arbitrarily prescribed, constant impedance. A complete trans- 
mission line formulation and Wiener-Hopf solutions have been 
described elsewhere 6). The geometry of the problem is illustrated 
in fig. 1. A plane wave is incident on the impedance half-plane at 
angle ¢; while the field is observed at a point P with the polar 
coordinates p, where —2/2 < ¢ < 3a/2. The angle at which an 
incident plane wave produces no reflected plane wave is denoted 
by da, the angle of ‘‘perfect absorption’’. All fields are independent 
of the y coordinate. Both kinds of polarizations are considered in 
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which either the total electric or total magnetic field may be along 
the y direction parallel to the edge of the half-plane. The rationalized 
MKS system of units and a suppressed time variation exp jw are 
employed. For each polarization the surface impedance Z of the 
absorbing half-plane is defined by the relations 


or (1) 
k 
Z = — COS dy 
if the total electric or magnetic field respectively is parallel to the 


edge of the half-plane. The propagation constant of the medium 
is k = w(ue)* where w is the angular frequency and mw and «é 


Tht 


z 
ABSORBING HALF -PLANE (x< 0,2z=0) 


INCIDENT 
PLANE 
WAVE 


Fig. 1. The geometry of the constant impedance half-plane. 


are the permeability and permittivity of the medium. In this 
way the impedance is normalized and specified in terms of a 
dimensionless parameter ¢q which has a simple and interesting 
physical significance. If the angle of incidence ¢; = ¢a, then a plane 
wave incident on an infinite plane surface with impedance Z is 
perfectly absorbed since no reflected wave of any kind is produced. 
If the same plane wave is incident on a half-plane with the same 
surface impedance, then it is “perfectly absorbed” in the sense that 
no reflected plane wave is produced, although a scattered, cylindrical 
wave is created. Williams 5) specifies the surface impedance in 
terms of an angle « which is equivalent to ¢g, but he does not 
introduce a physical interpretation of this angle as described here. 
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Our procedure allows us to specify the far-field scattering patterns 
in a physically significant manner for both types of polarization 
in terms of three similar, dimensionless parameters: the angle of 
incidence ¢;, the angle of observation ¢ and the angle of perfect 
absorption ¢yq. 

The Wiener-Hopf solutions provide the expected geometric optic 
fields (the incident and reflected plane waves) and a remaining 
diffracted field. This latter field may be evaluated by a standard 
saddle point evaluation to obtain the desired far-field behaviour 
(kp —> co). If the observation point P(p, 4) in the far-field is suffi- 
ciently far from the shadow-lit boundaries of the incident and 
reflected plane waves, then a pole singularity is sufficiently far from 
the saddle point. Under these simpler conditions the y component 
of the field for either polarization is given by 


F(p, $, bis ba) ~ Ald, $i, $a) f(p), 


where 
atkp 
and 
cos ¢; e?) ( = oT SS ba ) 
7 Ae sin 6 + sin da 
(P, Pir Pa) = 2(cos ¢¢ + cos da)(sin 8 — sin ¢4) 
t—- sin-@ | 
5 j S 
Ee cos 8 + cos da ( PASTY (3) 
] u2 u ] U4 Uu 
= ——_— a y —— oe d Vy 
(8) 2% Ju, Sin u ‘Ai 20 Jus. SiN u c 
uy = 6 + ba — %, 6=—¢ if —nz/2<$ <a/2 
Uz =04+ dg—G4, =az—¢ if a/2<¢ < 32/2, 
ug = 9 — da, e(z) = +.1 if z 2 0, 
us = bi — ha, 0 < (di or $a) < 2/2. 


The apparent discontinuity created by the e(z) term vanishes 
since cos § vanishes for z = 0. Moreover, to this order of approxi- 
mation the field on the imperfectly reflecting half-plane vanishes. 
The next higher order saddle point approximation is needed to 
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obtain non-vanishing (and discontinuous) fields on the imperfectly 
reflecting half-plane. Under certain conditions a pole singularity 
in the plane of integration at —z + ¢q will approach the saddle 
point if 6 - —4a and ¢q > 3a. This corresponds to a point of 
observation near the absorbing half-plane when the latter has a 
surface impedance which absorbs perfectly a plane wave at grazing 
incidence. Under these rather special conditions the results of (2) 
need to be modified. However, this correction is negligible for all 
values of surface impedance and angle of observation which we will 
consider. : 


§ 3. The diffraction patterns. To ascertain the properties of the 
far-fields the amplitude of the angular pattern function |A(¢, $i, ¢a)| 
given in (3) has been calculated as a function of ¢, the angle of 
observation. To cover a representative range of angles of plane wave 
incidence, the angles ¢; = 0°, 45° and 80° were chosen. For each 


-l00 -50 fe) 


50 100 
$-DEGREES 


Fig. 2. Diffraction patterns for different absorbing half-planes. 
(i = 0°) 


value of ¢; the representative angles of perfect absorption, ¢4 = 10°, 
45°, 80° and 90° were chosen. The resulting far-field diffraction 
patterns are shown in figs 2, 3 and 4. For each ¢, in each of these 
figures there are four curves, one for each of the four values of ¢q. 
One of these four curves (fq = 90°) is the well-known scattering 
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pattern for a perfectly reflecting half-plane when the magnetic field 
is polarized parallel to the edge of the half-plane. In addition, two 
more curves have been added to each figure. The first added curve 
(the dashed line) is the scattering pattern for a half-plane which 
is perfectly absorbing for ‘“‘angular propagating modes”; this is 
Sommerfeld’s black screen. Its far-field angular pattern function 
is given by 9) 
MLA 


m? — (x + bi — 4)? 


A(, $i) = : (4) 


ou 


.O1 


50 100 iSO 200 250 
¢ - DEGREES 


-100 -50 ° 


Fig. 3. Diffraction patterns for different absorbing half-planes. 
($i = 45°) 


The second added curve in each figure (the dotted line) is the 
scattering pattern for the well-known Kirchhoff black screen. 
For this type of black screen it is assumed that the tangential electric 
and magnetic fields vanish on the absorbing half-plane while the 
aperture field is assumed to be the incident plane wave. The re- 
sulting far-field angular pattern function is given by 1°) 


A(b, #) = —beot (S54). 5) 


The scattered fields contain the usual two shadow-lit boundaries 
for the incident and reflected plane waves at d = ¢iand¢ = a — $i 
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except when 4; = ¢q in figs 3 and 4. For these latter two curves the 
reflected plane wave is absent (“perfect absorption’’), and con- 
sequently the scattered field is finite and smoothly continuous 
across the non-existent shadow-lit boundary. As mentioned 
previously the scattered field amplitude tends to vanish as the 
observation angle ¢ approaches the surface of the imperfectly 
conducting half-plane. These factors serve to establish the general 
behaviour of the fields at the shadow-lit boundaries and the half- 
plane surface boundaries. 


10 


[A(e,4,.#,)| 


‘0.1 


oo $= 45° 
e e $= 10° 


——-—— SOMMERFELD BLACK SCREEN.| 6 
seers KIRCHHOFF BLACK SCREEN 


oo -50 10) 50 100 150 200 = 
 —DEGREES i 
Fig. 4. Diffraction patterns for different absorbing half-planes. 
. ($i = 80°) 


An interesting general feature of these curves is the presence of 
nulls in the scattering pattern for certain angles of observation as 
indicated in figs 2, 3 and 4. The angular pattern function A (¢, 4i, da) 
given by (3) vanishes if the terms in the square brackets vanish. 
For the range of values for d, ¢; and dq given in (3), these terms 
will vanish only if z < 0, that is if 7/2 < ¢ < 30/2. Consequently, 
nulls occur only in the back-scattered fields and not in the forward- 
scattered fields. With this restriction on ¢ the locus of the nulls is 
given by 

cos? 
eae (6) 

a 


For the special condition 4q = 4; when no reflected plane wave 


sind = 1 — 
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is present, (6) simplifies to ¢ = a — ¢;. The location of this zero 
coincides with the location of the pole singularity in A(d,4:, da) 
which gives rise to the geometric optic reflected plane wave. As a 
result, the reflected plane wave disappears for this special condition 
and the scattering pattern is a finite, smooth and continuous 
function of ¢ across the no longer existing shadow-lit boundary at 
ro) = 7% — di. 

When the magnetic field is polarized parallel to the edge of: the 
half-plane, then the half-plane is perfectly reflecting when ¢g = 2/2 
as indicated by (1). Under these conditions the well known scattered 
far-fields of the perfectly reflecting half-plane for all angles of plane 
wave incidence are symmetric with respect to the half-plane and 
exhibit a null in the aperture half-plane at ¢ = 2/2. For the 
absorbing half-plane, figs 2, 3 and 4 indicate that these symmetry 
properties are approached as ¢q approaches 7/2. Because of this 
property of the nulls, it is interesting and significant to interpret 
the asymmetric patterns of the absorbing half-plane as a kind 
of angular deformation or perturbation of the symmetric scattering 
patterns for the perfectly reflecting half-plane. This interpretation 
is valid for both kinds of polarization. Depending upon whether the 
electric or magnetic field is polarized parallel to the edge of the 
half-plane, we have a perturbation about a small admittance or a 
small impedance [¢q = 2/2 in (1)}. 

A study of the scattering patterns in figs 2, 3 and 4 shows that 
the scattered fields are most sensitive to the surface impedance 
in the back-scattering direction (z/2 < ¢ < 32/2) and much less 
sensitive in the forward-scattering direction (—2/2 < ¢ < 2/2). 
The half-plane with a surface impedance ¢g =’ ¢;, which produces 
no reflected plane wave, produces the smoothest scattering patterns. 
It is interesting to find that this “perfectly absorbing” half-plane 
produces approximately the same amount of scattering as the half- 
plane which is perfectly absorbing for “angularly propagating 
modes’’. This latter half-plane is Sommerfeld’s black screen. The 
calculations for this black screen obtained from (4) are much 
simpler than the corresponding calculations for the ‘constant 
impedance half-plane obtained from (3). In view of these results 
it may be plausible to suppose that the Sommerfeld black screen 
might serve as a simple, approximate mathematical model for a 
“perfectly absorbing”’ half-plane. 
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From the discussion given it should be easy to see how we might 
often obtain a reasonably accurate picture of the scattering pattern 
without actually performing the detailed calculations. To set the 
general level of the scattered field amplitude we might use the 
simple scattering pattern for the Sommerfeld black screen. Since 
the scattering near the shadow-lit boundaries tends to be independent 
of the surface impedance, the general behaviour near these regions 
could be reasonably approximated. Near the half-plane surface the 
vanishing of the fields could be estimated. Any nulls in the scatter- 
ing pattern could be located from (6). From all this information a 
scattering pattern function could be sketched which would exhibit 
the main properties of the diffracted field in a qualitative manner 
and which would not require elaborate calculations. 

I am indebted to Professors N. Marcuvitz and L. B. Felsen 
for their general guidance of this work and their many discussions 
and suggestions, and to Mrs. Toni Lotito for calculations of the 
scattered fields. 
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A “Wonderstone”’ microwave cavity for electron 
spin resonance spectrometers 


In our laboratory we use a 3-cm esr spectrometer with phase-sensitive 
detection (modulation frequency 1000 Hz). This spectrometer has been 
described in detail elsewhere 1). 

When using a brass TEo11 cavity we noticed that, even in the absence of 
electron spin resonance, the detected signal practically always differs from 
zero. This signal is caused by the interaction between the static magnetic 
field and the Foucault-currents induced in the walls of the cavity by the 
alternating magnetic field. This interaction causes a slight vibration of the 
cavity in a 1000 Hz-rythm so that the resonance frequency changes in the 
same rythm and the detected signal may contain a 1000 Hz component. 
This component is liable to fluctuate as a result of variations in temperature, 
klystron frequency and mechanical influences. These fluctuations manifest 
themselves as apparent noise, especially in measurements taking a long time 
as a result of broad resonance lines and application of long RC-times. 

To suppress this troublesome phenomenon, it is necessary to minimize 
the generation of the Foucault-currents, by making the conducting walls of 
the cavity as thin as possible Without weakening them mechanically. This 
can be achieved either in the way outlined by Eisinger and Feher 2), 
who, applying ENDOR technique, used a pyrex cavity coated with silver 
on the inside; or with the procedure suggested by Thompson et al 3) who, 
for the construction of their cavities with a low thermal expansion coefficient 
used silver-coated ceramic materials. Upon some orientations, however, we 
found that, at least for our laboratory, neither the coating of glass with a 
silver layer ensuring good adherence even at high temperatures, nor the 
machining of ceramic materials is an easy matter. 

A very easy solution however was obtained with the aid of the South- 
African ‘‘wonderstone”’ (mainly pyrophyllite), which had already been used 
before by other workers 4) e.g. as an insulating material. The properties of 
‘‘wonderstone’ have been described by Nel et al.5). If one follows the 
supplier’s instructions (Carters Ltd., London), the construction of a wonder- 
stone-cavity with the aid of conventional machine tools is not more difficult 
then that of a brass cavity. 

A section of a tunable wonderstone transmission TE 11 cavity, made by 
us, is shown in the figure. The outer section is rectangular, so that there is 
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sufficient material to screw the lids 1 and 2 and the brass waveguides 3 and 
4 on to the cylindrical casing 5. The cavity can be tuned by raising or lowering 
disc 7 with the aid of screw 6. Parts 2 and 6 are made of brass. 


CN 


iene ote 
aa 


Fig. 1. Section of a transmission TE 1; cavity. 1 and 2 lids, 3 and 4 brass 
waveguide, 5 cylindrical casing, 6 tuning screw, 7 disc. 


The inside of the cavity is coated with Dupont’s silver paste 6233. Without 
any polishing of the silver layer, the unloaded quality Q, of such a cavity 
amounts to about 12000. The Q, value of a brass cavity of the same design 
is about 10000, so that in this respect the wonderstone cavity is even better 
than a brass one. Moreover, the samples can be heated as easily ina wonder- 
stone cavity as in a brass-cavity 1). 

With the aid of the cavity described we suceeded in improving the ratio 
between the signal and the apparent noise by a factor 5. More detailed 
information about the construction of the cavity will be given elsewhere 8). 


Received 6th December, 1960. 
M. J. A. BAKKER 
J. SMiIptT 
Centraal Laboratorium 
Staatsmijnen in Limburg, 
Geleen, Netherlands 
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SHEAR FLOW OF A VISCOUS INFINITELY 

ELECTRICALLY CONDUCTING FLUID PAST 

A POROUS FLAT PLATE IN THE PRESENCE 
OF A TRANSVERSE MAGNETIC FIELD 


by A. S. GUPTA 


Department of Applied Mathematics, Indian Institute of Technology, Kharagpur, India 


Summary 


An analysis is made for a uniform shear flow of an incompressible, viscous 
and perfectly electrically conducting fluid past a porous flat plate in the 
presence of a transverse magnetic field. It is shown that in the case of 
suction an asymptotic profile similar to that of Griffith and Meredith (for 
the flow of an incompressible viscous fluid past a porous plate) exists if the 
suction velocity is greater than the Alfvén velocity. In the case of injection, 
however, the asymptotic profile will exist if the Alfvén velocity is greater 
than the velocity of injection. 


§ 1. Introduction. The flow of an incompressible viscous fluid 
past an infinite porous plate has been solved by Griffith and 
Meredith?!). Sakurai?) has solved the same problem when 
there is a uniform vorticity in the main flow. The present note 
is an extension of the problem considered by Sakurai to the case 
of a perfectly electrically conducting fluid in the presence of a 
transverse magnetic field. 


§ 2. Basic equations and theiy solutions. We take the x-axis 
along the plate (assumed to be non-conducting and non-magnetic) 
and the y-axis perpendicular to it. Since the plate is infinite, we 
can assume that all the physical variables are functions of y only. 

The equation of continuity is 

du/dy = 0, ) (1) 
which gives 
Vv = vo, (2) 


vp being the constant normal velocity at the plate. 
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The equations of momentum are 


du op d2u _ He alee 
== el: y 3 
ie dy Ox hice dy2 " 4% " dy . 
7 dH 
0 =— P ze Ay - (4) 
oy Ax dy 
When o = co, Maxwell’s equations give 
curl (q x H) = 0, (5) 
diy H = 0, (6) 
Since H, and Hy, are functions of y only, (6) gives 
HH, = constant — Ag, (7) 
and (5) gives 
H du a0 3 
—=0 : 
igre Cad (8) 
Since at y = co, u = U,, + Qy, we have from (8) 
HoQ = 9 
0 v0 dy (9) 
Equation (3) now gives 
op 
a ae 19) 
where 
eH? ) 
= p02 
a p ( Anpvy vo (1 1) 
Using (7), (8) and (10) we get from (3) 
du a du ow 0 . 
a? 2 dy et 
The solution of (12) subject to w = U,, + Qy at y = oo is 
u=U,+2 ciex beste l 
y+ Pp nO (13) 
Since u = 0 at y = 0, the above equation gives 
u=U, + Qy — U,, ex (-. 
y Nexen. (14) 
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In order that the above asymptotic solution may hold it is 
necessary that « > 0. Hence we must have 


Hel o* pu 15 

Avo es ( ) 
In the case of suction, v1 = — V, V > 0, the above inequality 
becomes 

V2 = V2 are, (16) 


where Vajz stands for the Alfvén velocity V (we/4p) Ho. In the 
case of injection, vo is positive, and the inequality can be written as 


a vo2, (17) 


which shows that the Alfén velocity must be greater than the 
injection velocity. 

In the absence of a magnetic field it has been shown by Sakurai 
that an asymptotic solution of the type (14) exists in the case of 
suction and not in the case of injection. However, the presence of 
a magnetic field in such a shear flow ensures the existence of 
asymptotic solutions of the type (14) both in the case of suction and 
injection provided the suction and injection velocity are respectively 
greater and less than the Alfvén velocity. 

From (8) and (14) we get 


H H U. 
d 2/04 BAe exp (— ay i 
dy vo wQ uQ 
Integrating this equation and noting that Hz = 0 at y = O (the 
plate being non-magnetic and non-conducting), we get 

Ho | | oy 

= — Dy + U,,| 1 — exp be I} (18) 
vo | 4 uQ 
Equation (10) gives on integration 


p = ax + Fly). (19) 


Ay 


Substituting (19) in (4) and integrating we get 
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Hence (19) gives the pressure distribution as 


PSO ed 02, (20) 
H, being given by (18). 
Received 21st January, 1961. 
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SOME PHENOMENA OCCURRING IN 
FLOATING BATTERY SYSTEMS 


by M. v. MENTS 


Research Council of Israel, Jerusalem, Israel. 


Summary 


The behaviour of a battery, its load and the supply system in floating 
condition is sometimes of a complicated nature. This applies in particular 
when the loading voltage is derived from a rectified a.c. supply which is not, 
or only partly, smoothed. In such a typical case it has been proved possible to 
express the equilibrium conditions of the system within a simple mathematical 
framework. The mathematical treatment and the results are elucidated and 
a numerical example met in practice is dealt with. , 


§ 1. Introduction. A battery continuously connected to a bus 
is said to float when the voltage of the charging line is slightly 
greater than the open-circuit voltage of the battery 1). There exists 
a great variety of examples where such a device is desirable. The 
conditions of the rectifier on the one side and of the load on the 
other put their stamp on the characteristics and the behaviour of 
the set-up. 

It is a known fact that battery rectifiers need not be smoothed 2), 
and the same is true for the rectifier in a floating system. In the 
latter case, however, some interesting phenomena occur which may 
give rise to more complicated situations and difficulties than in the 
case of simple loading. The complications arise in particular when 
the floating batteries have to deliver current during a part of the 
cycle. 

Such a case was encountered in realising the tension supply of the 
National Standard Clock in Israel 3). Because of the frequently 
occurring mains supply failures, use had to be made of a floating 
battery system for both high and low voltages. It transpired that 
equilibrium conditions were obtained after 5 or 10 hours only, and 
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that there are situations where the floating voltage at the end is 
either lower (which means that the battery is slowly discharging) 
or higher than the prescribed values *). 

In order to be able to handle the situation and to know in advance 
what values should be used for the several resistances and voltages 
of the circuit, a simple mathematical treatment of such a system 
was developed. The results proved to be very successful in promoting 
a better understanding of what happens in such systems. Some 
approximations and assumptions had to be made, so as not to 
complicate the calculations more than necessary. It has been 
verified, however, by practical observations that this mathematical 
treatment gives very reliable results. 


§ 2. The general set-up of a floating battery system. A general 
scheme of a floating battery system is sketched in fig. 1. The 
internal resistance of the rectifier is given by R,1 and, as will be 
explained, may have to be completed by an extra one R,2. The 
internal battery resistance Rp; will be greater by an amount Koz, 
as caused by relay-contacts, wiring, etc. +). The load resistance Ry 
may be greater by the amount R72 for the same reasons, or by 
purposely adding resistance. 


Rr Rez 
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Fig. 1. General set-up of a floating battery system. 


For a given time dependence of V,, all these quantities are inter- 
related in some way, and some of them determine the value of the 
others. In most cases, part of the quantities may be fixed and 
invariable, right from the beginning (e.g. V;, Ryi, Ri), another part 
may have a prescribed value (say: the average load current é3 and 
the battery voltage V» in floating position) and the rest (e.g. Ry) 


*) See footnote in § 3. 
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may be chosen within certain limits. In the next section, the re- 
lations between the several quantities will be evaluated. 


§ 3. Mathematical treatment. Let 


Rr ay Rye tr Ry, 
Roi + Roz = Ro, (1) 
Ra + Riz = Ry 

and further 


R,/Ri =; (2) 
Ry/Ri = 6, ©) 
RrRi + RpRi + RrRo = 7h bot yb e\s (4) 


R? 


If it is ignored for the moment that 7; cannot be negative in the 
case of an ideal rectifier, then the Kirchhoff relations 


11 = 12 + 123, (5) 
Ve = 11R,y —= 1oRp =— Vo = ‘O} (6) 
igsR, — tgp — Vp = 0 (7) 


should yield the currents in the several circuit branches. The result 
is found to be 


(1 + 0) V; — Vo 

— 8 
A sR : ( ) 
4 Ve Ni) Ve 

— 9 
12 sR, , ( ) 
: bV~+1V>» 

= SO. 10 
13 Rs ( ) 


However, for an ideal rectifier, 7; cannot be negative, so that 
(8), (9) and (10) only hold when V; > Ve, where from (8) the critical 
voltage Ver is given by 


Vo 
= ———_-,, l 1 
Ver= 35 (11) 
During that part of the cycle for which V; < V er there holds 


Vo 
Ri(l + 8) H2 


44, = 0, te = 13) = 
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It will further be assumed that use is made of a full wave rectifier 
and that V; is given by 
V, = |V, sin wt}. (13) 


When no inductive or capacitive elements are present in the 
whole system, only the phenomena taking place during the time 


OS Cape? eT Foeseteal 


Fig. 2. Full wave rectification. 


periods 0 +¢, and-¢; ~ ¢g will have to be studied (see fig. (2)), 
where 


tga 4] = wf 2e. (14) 
The time ¢; for which V; = Ve, is defined by 
Ver = Vr sin ott, (15) 
and using (11) it follows that 
sin wt, = a dt Me . (16) 
(1 + )P, 


It will prove to be useful to define an auxiliary angle « by 
a = 2/2 — oty. (17) 


Now it is assumed that the battery delivers the load current during 
that part of the cycle for which V; < V,, and behaves then as if 
no rectifying system were present at all. As a result, the battery 


loses (between ¢ = 0 and t = ¢;) an amount of charge given by 
th 


You fv a Ran b) (5 se *); 8) 


0 
as can be easily verified with the help of (12) and (17). 
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The amount of charge, received by the battery during the second 
part of the cycle (between ¢; and fz) is given by 


te m/2 


Qin = | 72 di ~| 


ti a/2—« 


V, sin wf — Vo(1 + 7) 


dot. 
osk 


The result of the integration is 


V, sin x Vol + ro 


wosRy wosk, 


Qin . 


(19) 


It is assumed that Ry; has the same and moreover constant value 
for both the time intervals 0 — ¢; and ¢; — ¢g. This is true only in 
first approximation, but as in any case Ry <R; in practically 
occurring systems, one need not bother about deviations. This will 
find its expression in (24). 

The net charge per cycle is given by 


Qe har 4(Qin — Qout) 
and the net charge per second by 


] : 
Qs = - Qe = te. (20) 
Combination of (14), (18), (19) and (20) yields 
om(1 + 6) Rig _ sina(1 + 8) eo Ul +r)(l+d)e 2 


Vo a sVo s 2 Azle 
By making use of (2)—(4), showing that 
(I+ 6)(1 +7) =s +1, 
and of (16) and (17), yielding 
V(l+6) 1 
Vo cos a” 
it is easily seen that (21) goes over into 
jas(1 + 2) = tga —« (22) 
where ¢ is defined by 
a) Bee (23) 


Vo 
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The value of ¢ is directly connected with the practical considerations 
that prescribe the value of V» *). 

Although no well-defined data can be found in literature with 
respect to the residue charging current 72, it has been proved by 
studying the instructions of the manufacturers that the charging 
current in floating condition is very small as compared with normal 
charge or load currents *). The electro-chemical processes occurring 
in this situation, where the battery tries to pull itself up to reach 
asymptotically the average bus voltage, play an important part 
of a very complicated nature. 

In normal situations where b <1 it follows from (23) that é is 
given by € & 79/73, so 

e <= I. (24) 


(Only in cases where for some reasons 6 > 1 it follows that 
é w 12Rp/Vo, 


and it might be that the value of ¢ is then comparable with or 
even greater than unity). Inserting (24) into (22), yields 


tas ~ tga — a. (25) 


It seems that (25) (or eventually (23)) governs the relationship 
between all the quantities appearing in floating battery systems. 
This result will be demonstrated with an example, met during the 
installation of the National Standard Clock in Israel 3). Here the 
following quantities were encountered for a low voltage supply: 
Vy = 6.5 V (three cells; value given by manufacturer). 

Ry = (0 OF Rie ==(0S Q, Ru = 2.0 0 (average, if Vu = oR | Va 

43° =O0:107A; 

and use was made of a full wave selenium rectifier for which 

V, = 15.7 V, Ry = 0.70 Q, as measured with dummy loads. 
Now the following values can easily be evaluated: 


b == 0.093 &. ==-0.036 
sin wt, = 0.38 So a DTG 

a == 11 6radian yp =n 0.61 
tg « — a = 1.26 Ry i Yoon 


*) In practice, Vy varies between 2.15 and 2.20 volt, depending on the type of battery; 


the current 7g is of the order of 0.5 to 1.5% of the 8 hour rate of the battery. 
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As Ry1 is 0.70 Q, it easily follows from (1) that an extra resistance 
Rrg = 0.61 Q has to be added. 


§ 4. Conclusions. With the help of the framework previously 
developed one can define the mutual relationship between the 
various parameters of the system and calculate their values. 

It further enables us to investigate what happens in the case of 
slight variations of one or more of the values with respect to equi- 
librium values. Thus it is easily seen that 


Vo nV, ( cos2a ) 


(26) 


—— L p)2 
OR, BRpeelis 
which means that one immediately knows how to change the added 
resistance R,2 in case of a drop or a rise in the observed battery 
voltage as measured at the terminals. Such events might occur e.g. 
by the ageing of a selenium rectifier. (In fact, it was experienced 
in the case of the example dealt with in the previous paragraph 
that the voltage dropped by 0.20 V after a couple of months. The 
situation could easily be restored by reducing the resistance Ry 
by 0.10 Q as follows from (26). A new equilibrium was obtained after 
a couple of hours, Vj having reached again the prescribed value of 
6.5 V.) 

Another practical application is to see whether the nature of 
(23) still offers enough freedom to cope with one or more other 
conditions. In the example treated it was e.g. highly desirable that 
the average load current 73 should remain the same in the case of 
disconnecting the battery. It can be shown that this condition 
implies 


sin « 


1 I be 

=u 4a(i+ 5): ro 

With the help of (22), (23) and (27) one is able to solve for all the 

quantities, to find amongst other things that V, should have been 
chosen 30 V instead of 15.7 V, and R; = 4.1 Q. 

The average load current would have been the same (3.05 A) 

with disconnected battery, whereas in the present case this current 


has a value of 


anf, 2 Vy 
emery) 
if one disconnects the battery. 


= ILE 
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Finally it is worthwhile investigating the ripple on the load 
voltage. Fig. 3 is drawn with the help of the results of (10), (12) and 
(13) and the values of the numerical example. The real time de- 


wl, TI TT 2 oda en 


Fig. 3. The ripple on the load voltage; the figure is drawn according to the 
values of the numerical example. 


pendence could be verified by making use of an oscillograph and by 
measuring the voltage drops over the several resistances. The 
agreement between calculated and observed data provides a 
justification of the method as a whole. 


Received 28th January, 1961. 
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AXISYMMETRIC STAGNATION POINT FLOW IN 
MAGNETOHYDRODYNAMICS 


by W. Il. AXFORD 


Defence Research Telecommunications Establishment Ottawa, Canada *) 


Summary 


This paper contains a theoretical treatment of the axisymmetric stagnation 
point problem in magnetohydrodynamics. Two distinct arrangements of the 
magnetic field are considered, with the lines of force lying in the plane of, and 
perpendicular to the streamlines, respectively. The boundary conditions 
are examined carefully and the solutions are completed by joining them to 
meaningful distributions of magnetic field in the solid region. Exact solutions 
can be obtained by numerical methods, but it is more illuminating to consider 
approximations valid for very large and very small diffusion numbers (ratio 
of kinematic viscosity to magnetic diffusivity). Distinct boundary layers of 
current and vorticity occur in these limiting cases, which can be treated 
separately. For flows with the first arrangement of the magnetic field, and 
in which the field does not vanish at infinity, it is deduced that the solution 
refers to a forward stagnation point when the Alfvén number is less than 
unity, and to a rear stagnation point when the Alfvén number exceeds unity. 
An interesting choking effect is found in flows with the second arrangement 
of the magnetic field when the diffusion number is small, and the Alfvén 
number exceeds unity; it is believed that this effect could be demonstrated 
in the laboratory. 


§ 1. Introduction. A number of papers have been written in the 
past few years on the subject of magnetohydrodynamic boundary 
layers. Most attention has been paid to the case of weakly conducting 
fluids, with the aim of predicting the effects that may be found in 
aeronautical situations, although other limiting cases have also 
been considered. In the present paper we hope to clarify some of 
the more general aspects of the subject without emphasizing particu- 
lar numerical results. The problem we shall consider is that of 


*) Work carried out at Cambridge, England. The author is a member of the New 
Zealand Scientific Defence Corps. 
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axi-symmetric stagnation point flow of an incompressible, viscous, 
electrically conducting fluid in the presence of a magnetic field; 
this problem has the advantage that exact (numerical) solutions 
can be found if necessary to compare with approximate solutions. 
The term ‘boundary layer’ is used in a relative sense in this context 
and does not presuppose the existence of large Reynolds numbers. 
Several authors have considered the stagnation point problem in 
magnetohydrodynamics, both for the two-dimensional and axi- 
symmetric cases, the first being Michael!) who examined the 
current layer in the two dimensional case with viscosity neglected. 
Michael’s solution demonstrates the manner in which joint con- 
vection and diffusion of the lines of force produces a boundary 
layer of current, but as the existence of similarity solutions was 
not recognized by Michael, mathematical difficulties forced him 
to consider only the case of vanishingly weak magnetic field. 
Similarity solutions of the two-dimensional problem were examined 
by Neuringer and McIlroy 2), whose methods, however, were 
questionable. This work was criticized by Hess 8), but his paper 
did not succeed in reducing the confusion to any great extent. 
The impression was given in Neuringer and MclIlroy’s work 
that their solution is exact, leading Rossow *) to believe that the 
close agreement with his approximate solution confirmed the 
validity of the approximations. However, as pointed out by Hess, 
the Rossow approximations were made by Neuringer and 
McIlroy in a hidden manner and the agreement of the solutions 
was to be expected. Rossow’s boundary layer equation ®) should 
be used with caution, for although it gives correct results in some 
cases 6), there are other cases in which a careless treatment leads to 
results that are meaningless in the sense that they cannot be related 
directly to any flow problem. This is exemplified in the papers by 
Neuringer and McIlroy and Rossow 4), where the behaviour 
of the magnetic field and the boundary conditions it must satisfy 
are not considered in any detail; it is not surprising then that we 
find that the results obtained by these writers are not particularly 
relevant to aeronautical situations (Hess). It is true that for 
a weakly conducting fluid, induced fields can be neglected in the 
viscous boundary layer, but this does not imply that induced fields 
can be neglected entirely as the field may interact strongly with 
the flow elsewhere, causing the field in the boundary layer (and 
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within the boundary as well) to be distorted so that it cannot be 
specified in an @ priori manner. The correct method of attacking 
the problem is given in this paper; it has also been given by M eyer’), 
who has obtained some numerical results. Further numerical results 
have been given by Poots and Sowerby 9), 

In cases where a shockwave occurs a short distance upstream of 
the stagnation point, and the fluid is electrically conducting only 
between the shock and the surface of the body, induced fields can 
be neglected if the magnetic Reynolds number based on the stand- 
off distance is small. Several authors have considered problems of 
this type, notably Meyer 8), Kakutani 9) and Oguchi 1%), 

In contrast to the papers referred to so far, which treat only the 
case of weakly conducting fluids, the present paper treats both 
limits of very large and very small electrical conductivity. Two 
distinct arrangements of the magnetic field are considered; in the 
first arrangement the lines of force lie in planes which intersect the 
axis of symmetry, while in the second arrangement (which has not 
been considered before) the lines of force are circles and result from 
an applied axial current. The Alfvén (S) and diffusion *) (¢) numbers 
are found to be the most important parameters in the problem. 
For small diffusion numbers the region of flow in which diffusion 
of the lines of force (current layers) occurs is wider than the viscous 
boundary layer by a factor O(e—?), whereas for large diffusion 
numbers the viscous boundary layer is wider than the current layer 
by a factor O(e+). An interesting effect is found to occur in flows 
with the first arrangement of the magnetic field when the diffusion 
number is very large; whereas solutions for only the forward 
stagnation point are possible in ordinary hydrodynamics, the 
solution in the present example refers to a forward stagnation point 
when the Alfvén number is less than unity, and to a rear stagnation 
point when the Alfvén number exceeds unity. It is believed that 
the same effect occurs when the diffusion number is small, but 
confirmation from numerical solutions of the equations is required. 
A similar effect was found by Carrier and Greenspan }%) in their 
_ *) The term “diffusion number” is used throughout this paper to denote the ratio of 
kinematic viscosity to magnetic diffusivity. It seens to have been Batchelor 11) who 
first skowed the significance of this quantity in magnetohydrodynamics. Sooner or later 
the quantity (magnetic diffusity/thermal diffusivity) must occur in the literature, 


and so, to avoid confusion it is suggested that the term “magnetic Prandtl number’’ 
should not be used. 
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analysis of the flat plate boundary layer problem. In flows with the 
second arrangement of the magnetic field, ‘choking’ occurs when 
the Alfvén number approaches unity and the diffusion number is 
very small; it is suggested that a laboratory demonstration of this 
result may be possible. 

It should be noted that the technique of separating the current 
layer and the viscous boundary layer, when the diffusion number 
is both very large and very small, cannot in general be regarded as 
forming a first approximation for a series solution in ascending or 
descending powers of ¢, although one further iteration may be 
possible. However solutions can always be obtained in the form of a 
series in ascending powers of S; when S = 0, the solution for the 
velocity field is that obtained by Homann !}8), and the structure 
of the current layer is similar to that of a thermal boundary layer 
in the absence of buoyancy forces }4). 

Previously, many writers have ignored the details of the magnetic 
field in the solid regions, usually stating that a distribution of 
surface currents can always be found so that the conditions at 
the boundary are satisfied. In the author’s opinion, a problem of 
this type should always be completed by joining the solution in the 
fluid to a suitable solution for the solid region 15), otherwise the 
feasibility of the problem, as well as its relevance to any physical 
situation, may be in doubt. The solutions given here do not depend 
upon whether the solid is electrically conducting or not, and in 
particular they represent the steady state in all cases for which the 
conductivity of the solid is finite. Some authors (notably Michael) 
have made the assumption that the solid is a perfect conductor, 
but this is hardly reasonable for real materials in steady flow 
problems unless a superconductor is involved. For example, a rough 
calculation shows that any distortion of the magnetic field at the 
boundary of a solid copper body penetrates a distance of 10 centi- 
metres in a time of about 2 seconds, which is not long in engineering 
contexts, except perhaps in shock tubes. 


§ 2. Stagnation point flow with no azimuthal magnetic field. When 
there is axial symmetry, it is convenient to use cylindrical polar 
coordinates (7, 0, z) such that the stagnation point is at the origin 
and the regions z > 0 and z < 0 are occupied by fluid and solid 
material of uniform composition, respectively. The equations of 
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magnetohydrodynamics in these coordinates have been given by 
Yih 16) and others; for brevity we simply note that the solution 
of these equations for stagnation point flow, when the velocity and 
magnetic induction vectors have the forms (w, O, w) and (B;, O, Bz) 
and there is axial symmetry (0/20 = 0), is 


yp = Ur/(«/U) {($), A = Br? y/(«/U) g(¢). (2.1) 


Here y is a stream function and A the non-zero component of the 
magnetic vector potential, such that 
1 op 1 dp 1 eA. 1 0A 


“= — — , w = — ——, By = Ne Se 
bet te Mees has Byte BR : i a 


the independent variable is defined to be 
b= V(Uln):, (2.3) 


and the constants U and B have the dimensions (velocity/length) 
and (magnetic induction/length) respectively, while « has the 
dimensions of a diffusivity. 

The pressure # has the form 


b = po — 3pU**hi(p) — 2pUKpa(4), (2.4) 


where p is the density of the fluid and fo is a constant. The functions 
(4), g(¢), p1(¢) and p2(¢) satisfy the following ordinary differential 
equations: 


~ f” + 2ff — f? + pi = See", (2.5) 
alg” = 2g} — re", (2.6) 

pr’ = S(e’)', (2.7) 

Ps = (7 +—1", (2.8) 


where primes denote differentiations with respect to ¢, y is the 
kinematic viscosity and 7 the magnetic diffusivity of the fluid, 
S = B2/4nypU? is the Alfvén number, and yw is the permeability 
of the fluid. The electric current and the vorticity have only 
azimuthal components, being 


Dee ~ V/(U|«) 2°(9), 29 = Urv/(U]x) (9) (2.9) 
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respectively; the electric field vanishes everywhere. As in the 
ordinary hydrodynamical problem (the equations for which are 
found by taking B = S = 0 and« = y in the above), the equations 
determining the pressure can be integrated at once; in particular 


Ppa yg eh yl, Sy (2.10) 


where we have written the constant of integration as (1 — Sa). 
Substituting this result, (2.5) becomes 


~~ #" + 2if” — f+ 1 = Sl2ge” — g* + 2). (2.11) 


We are at liberty to choose x in any way we please, but as we are 
interested only in approximate solutions of the equations for very 
large and very small values of the diffusion number, it is convenient 
to take xe = v when considering layers in which viscosity dominates 
(i.e., viscous boundary layers), and (usually) « = 7 when considering 
layers in which diffusion of the lines of force is important (i.e., 
current layers). For the viscous layers we use F, G and & in place 
of f, g and ¢, and the equations to be solved are therefore 


PIF ROL een 0G es See (2.12) 
G" — 7G — BG): (2.13) 


similarly, in the current layers where we use ¥, Y, and #@ in place 
of /, g, and ¢, the equations become 


ef’ ae LFF 2 sa F'2 — j= S(2g 9 = G'2 + a), . (2.14) 
G" — DF'G —~ FG". (2.15) 


Note that dG/d§ = d¢Y/d@ and G(é) = e*Y(9), with similar 
relationships for F(&) and #(#). For numerical integration of the 
full equations, it may be convenient to choose « = 4/(v7). 


Boundary conditions. For stagnation point flow, the boun- 
dary conditions. to be satisfied by the velocity components are 
“u=v=0O0onz=0. and u > Ur as z > oo; thus 


{(0) = f'(0) =0, and f’($) +1 as 6 co, (2.16) 


and we see from (2.6) and (2.9) that the electric current vanishes 
on the boundary .The tangential components of the magnetic 
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induction must be continuous at the surface z = 0. Hence, if dashes 
are used to denote quantities in the solid region z <0, then onz = 0 


B’,|u' = Brg'(0)|u, B’z = —2B4/(x/U) g(0), (2.17) 


and the electric field is zero. We shall consider the case in which 
no currents flow in the solid region, so that the magnetic field 
is curl-free; a general solution which enables the boundary con- 
ditions to be satisfied is 


A’ = (az + 0) Pr, (2.18) 


where a = Bg’(0)u'/w and b6 = By/(«/U) g(0). The only solution 
for which there is no field in the solid and the conductivities are 
finite is the trivial case of zero field everywhere. There are two 
possible ways in which the conditions on g(¢) can be chosen; 
these are: 


Case A g'(¢) > 1 as d > oo, g(0) = Dv/(U/k), (219) 
Case B g(p) +0 as d +00, g(0) = Dy/(U Jk), (2.20) 


where D is a constant depending on the form of the magnetic field 
as z — — oo. These conditions determine the value of «, for by 
substituting the values at infinity in (2.11) we find that « = 1 
in case A, and « = Oincase B. 

The stagnation point solution gives a locally valid representation 
of a large number of possible flows. The solution in case A could 
represent, for example, the flow near the stagnation point of a 
sphere which is moving at constant speed along the lines of force 
of a uniform ambient magnetic field. Similarly the solution in case 
B could represent the flow near the stagnation point when a 
magnetised sphere moves at constant speed in the direction of its 
magnetisation, the magnetic field vanishing at large distances from 
the sphere. In either case the sphere may contain, say, a suitably 
aligned dipole at its centre, and (2.18) remains a locally valid 
solution for the field near the stagnation point. The only condition 
that must be satisfied is that the boundary is effectively plane at the 
stagnation point; this is true provided both »/(y/U) and +/(»/U) 
are small compared with the radius of curvature of the body at 
the stagnation point. 


Approximate solutions valid for small diffusion 
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numbers. For very small values of the diffusion number, viscous 
effects are negligible in the current layer, so that the appropriate 
equations for this layer are (2.14) and (2.15) with the term eF” 
omitted from (2.14), thus 


LFF” — F'%*4+ 1 = S(2GF" — G2 + a), (2.21) 
and 
G" — DF'G— DFG’. (222) 


A relatively thin viscous boundary layer occurs next to the boundary, 
permitting the no-slip condition to be satisfied; this condition is 
ignored when the above equations for the current layer are solved. 

It is not easy to obtain a solution of (2.21) and (2.22) in any given 
case owing to their behaviour at ¥ (i) = 0; however, approximate 
solutions can be sought in series form: 


GO oot. F (P) = > SF ale). 
0 0 

Clearly Fo(%) = # in all cases, and Yo(*#) satisfies the equation 
Go" + 20G'y — 299 = 0, (2.23) 

which has the solution 

B co 
Go(9) = «oP + G(0)(1 — 2ffe™ dx dy). (2.24) 
Oy 


Since Y’9(0) = « — a#GY(0), the solution in the solid region to this 
order of approximation is given by 


a= 2 | («= DY) e+ wm’ | +019) (2.25) 


The reader will note that the field in the solid region can be uniform 
(as assumed by Neuringer and McIlroy) only in case A when 
D = y/(n/xU) if S= 0. Solutions for the higher order terms are 
best obtained numerically if required, but it is interesting to 
examine the equation for #1(#) which becomes on differentiating 

FF") = —2GoG 0"; (2.26) 


this can be integrated to give 


(0) = MJ Golx) e~™ ax, (2.27) 
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where M is a constant. We see that #”j(0) is finite here, although 
the vorticity is often logarithmically infinite in the O(S) solution 
for this type of problem (see for example, Ludford andMurray =) ye 
The solution for Yo(#) with « = 1 and D = 0 has been given pre- 
viously by Meyer) and Poots and Sowerby 19); (2.24) and (2.25) 
represent the general solution for arbitrary D and « = Oor 1. 

Compared with the current layer, the viscous boundary layer is 
so thin when e <1, that as a first approximation the components 
of the magnetic field can be taken to be constant; however a certain 
amount of care is required in order to obtain the correct approxi- 
mation for the Lorentz force terms in (2.12). From (2.15) we note 
that in the viscous boundary layer 


Oar ‘ d2gG dF 


= a 2e*G (0) = De F'(0)-F, 2.28 
ag ES Os 2690) (2.28) 
and so 

dg dF 
ale) ey 2etG(0) —. (2.29) 
dé \ dé dé 


Thus the change in dY/d# across the layer is O(e?), confirming our 
assertion that the magnetic field is effectively constant in the viscous 
boundary layer. Hence we see that, correct to O(e?), 


2GG” — G2 = 4G(0)2F’ — F'(0)?, (2.30) 


and so an approximate equation for the viscous boundary layer, 
valid when ¢ <1, is 


FY” + 2FF” — F’24 1 = S(4G(0)2F’ — G'(0)2 + a). (2.31) 


The boundary conditions are F(0) = F’(0) = 0 and F’(é) > E as 
£ — oo, where E is the positive root of the quadratic equation 


E2 + 4SY(0)2E + S(« — Y’(0)2) —1=0. (2.32) 


The equation used by Rossow 4) is a special case of (2.31) with 
G'(0) = 0 and « = 1 (that is, case A), although it is by no means 
clear from Rossow’s presentation that this is the problem being 
considered. Hess 3) suggested that the magnetic field in this case 
would have to arise from a magnet placed in the stream ahead of 
the stagnation point; this is not a correct interpretation of case As 
but it does emphasize that this case is not particularly relevant to 
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aeronautical problems such as re-entry. Clearly there is a need for 
careful approximations in these problems, since it is doubtful 
whether the reasoning employed by Rossow 5) would lead to (2.31) 
for general values of Y(0) and Y'(0) with « = 0 or I. 

In order to solve (2.31), it is first necessary to solve (2.21) and 
(2.22), so that a value can be obtained for Y’(0). Normally one 
solves equations of this type numerically by integrating towards 
the origin (A =), starting with suitable asymptotic solutions 
valid for large values of #. In this case, the procedure is complicated 
by the tem FF" in (2.21), which contains the highest derivative 
of ¥, but which vanishes at # = 0. A second series solution, valid 
for small values of 3, must therefore be used in place of the numerical 
solutions when # is small. (Meyer “) has obtained some numerical 
solutions for case B, but does not make this point explicitly). 
This procedure is not always possible in boundary layer problems 
of the type U = cx™, as the equations for the current layer have 
a non-integrable singularity at the boundary when m < 0 (Lewel- 
len 2°)), With Y'(0) known, the solution of (2.31) is straightforward, 
although one must make use of an asymptotic solution valid for 
large € to start the integration. 

This approximate method of solving the stagnation point 
problem for small values of ¢ is considerably easier than attempting 
to solve the full equations (2.6) and (2.11) as they stand. The 
difficulties encountered in this problem are largely caused by the 
two-point boundary conditions. If, for example, one tries to solve 
the full equations by integrating towards the origin, using asymptotic 
solutions to start, the correct solution must be chosen from a two- 
parameter family of solutions since the conditions at the boundary 
leave two arbitrary constants in the asymptotic solutions. The 
approximate method avoids the difficulty of having to estimate 
these two parameters concurrently by treating the problem in two 
steps (that is, as two boundary layers) with one parameter to be 
determined in each. 

For very large values of S, an approximate solution of (2.31) can 
be found by neglecting the terms other than F” on the left hand 
side. The resulting equation is linear and its solution is easily found 
to be 

G'(0)2 — « 


Ree 4a mr | — exp(—25'9(0) £1. (2.33) 
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It can be shown that this is a Hartmann type of solution 21), with 
the boundary layer thickness inversely proportional to St. Note 
that the solution indicates that reversed flow can occur in case 
A if S > 1; this effect is referred to again at the end of the section. 
In case B it is likely that no solution exists when S > 1, since the 
appropriate solution of (2.21) with the left hand side neglected, 
cannot be made to satisfy (2.20). 


Approximate solutions valid for large diffusion 
numbers. When the diffusion number is very large, we assume 
that the flow in the viscous boundary layer is, to a first approxi- 
mation, that of a perfectly conducting fluid, and the adjustment 
of the magnetic field to satisfy the boundary conditions at z = 0 
takes place in a narrow current layer next to the surface. The 
current layer is not thinner than the viscous boundary layer by a 
factor e~? as might be expected to begin with, but a certain negative 
power of ¢ is involved, and our assumption that the current layer 
is negligibly thin for sufficiently large e can be justified a fosteriort. 

For the viscous boundary layer we neglect the term G”/e in (2.13), 
which can then be integrated at once to give 


== ol (2.34) 
This result can be substituted into (2.12) to give a single equation 


for F(&); in case B this is simply the ordinary Homann equation, 
and in case A the equation is 


F” + (1 — S)(2FF” — F’2 + 1) =0. (2.30) 
If S < 1, we can make the transformation 
F(Ey= (1 —"S) #68), = (1 1S)i6, (2.36) 
so that (2.35) also reduces to the Homan equation 
0” + 200” — 6%+1=0, (2.37) 


with boundary conditions 0(0) = 0’(0) = 0, (co) = 1. Note that 
the transformation shows that the effect of the magnetic field is 
formally equivalent to an alteration of the kinematic viscosity from 
» to v/(1 — S) *) if S < 1; the viscous boundary layer thickness is 
correspondingly greater by a factor (1 — S)*+. 


*) F. A. Goldsworthy (private communication). 
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When the Alfvén number exceeds unity, there is no solution of 
(2.35) if the flow is directed towards the boundary *), but there 
is a solution if the flow is reversed. That is, taking U to be positive 
as before if w + — Ur as z > 0, (2.16) is amended to /’(¢) > —1 
as # > oo. We regain the Homann equation (2.37) once again by 
making the transformation 


F(é) = —(S — 130), y= (5 — 1 )ie, (2.38) 


in (2.35); the boundary conditions, as before, are O(0) = O’(0) = 
= 0, O’(co) = |. In this case the effect of the magnetic field is 
equivalent to altering the kinematic viscosity to »/(S — 1), thus 
changing the thickness of the viscous boundary layer by a factor 
(S — 1)~}. There appears to be no solution of the stagnation point 
problem when S = 1. 

It is assumed that the current layer is so much thinner than the 
viscous boundary layer that the velocity profile in the current layer 
is effectively linear; thus (2.15) can be written 


G’ + 28929’ — 480g = 0, (2.39) 


if F (3) ~ Pd in this region. In case B, it is found that 6 = 0.656e7}, 
and in case A 


B = 0.656(1—S)t/e#, (S <1); B =—0.656(S—1)#/e!, (S>1). (2.40) 


To obtain the correct scale for the current layer we must make a 
change of independent variable to 


A = Bd, (2.41) 
so that (2.39) becomes 
ae + 2A2 nah 41g =0 
2 Ter ==/0) (2.42) 


Thus the scale of the current layer is not 4/(7/U) as it is when 
é <1, but instead it is 6-4y/(y/U). The ratio of the thickness of the 
current layer to that of the viscous boundary layer is O(e—*), which 
justifies our original assumption that this ratio is small for large 
diffusion numbers. Correspondingly, the ratio of the term neglected 
in (2.13) to those retained is G’/eGF’ ~ O(1/e&8) near € = 0, which 


*) Since the Homann equation cannot be solved with the boundary condition 


O’(co) = —1. 


STAGNATION POINT FLOW IN MAGNETOHYDRODYNAMICS 225 


is not small when & = O(e-*) and indicates the necessity for a 
current sheet at the boundary. By substituting « = y/8* in (2.11) 
we see that d°F /d/3 = O(e-122), and hence the solution F(9) = p92 
is valid in the current layer provided 4 < e!. 

Note that (2.42) has the asymptotic solution 


Pinay fart? datos 01), (2.43) 
Aa 


for large 2, where Ny and N¢2 are constants. Clearly Nz = 0 in case 
B since Y(A) +0 as 4 +00; in case A, ¥~ F as A > 00, So 
that No = Bs. 

The behaviour of the solution for case A is particularly interesting. 
It is a feature of magnetohydrodynamics that vorticity is not only 
diffused by viscosity but is also propagated along the lines of force 
as magnetohydrodynamic waves. In this example the vorticity 
produced as a result of the no-slip condition is diffused by viscosity 
to form a boundary layer and at the same time tends to be propagated 
away from the surface. When S < 1, the speed of propagation of 
magnetohydrodynamic waves is less than the fluid speed, and so 
the vorticity remains confined to the neighbourhood of the bounda- 
ry, although the tendency to propagate is made evident by the 
increase in the boundary layer thickness. When S > 1, the vorticity 
can move freely against the stream and so does not remain in the 
form of a boundary layer when the flow is directed towards the 
boundary: in other words the boundary layer does not exist 
because there is no longer any balance between convection and 
diffusion of vorticity. If the stream direction is reversed, however, 
the vorticity which moves against the stream remains next to the 
boundary, and the boundary layer is kept intact; in this case the 
solution refers to a reversed stagnation point, and the flow away 
from the boundary represents a wake. 

The vorticity argument given here does not apply when the 
diffusion number is very small, and so it is not immediately clear 
whether or not the reversed stagnation point flow is possible in 
this case. Solutions of (2.21) and (2.22) appear to be possible when 
the condition at infinity is amended to ¥’(co) = —1 and both 
roots of (2.32) are negative in case A if S > 1/[1 — 9’/0)?] > 1, 
so that a solution of the viscous boundary layer equation (2.31) 
should also be possible for reversed flow. A final answer to this 
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question must be left until an attempt is made to solve the equations 


numerically. 


§ 3. Stagnation point flow with a purely azimuthal magnetic field. 
When the magnetic field is produced by an axial current and is 
therefore of the form (0, Bo, 0), a solution of the stagnation point 
flow problem can be found with the stream function p defined as 
in (2.1) and with the magnetic induction given by 


Bo = Brh(d), (3.1) 


where B is a constant with dimensions (magnetic induction/length). 
The total magnetic and fluid pressure vy has the form 


4 =p + Bo? [Buu = xo — ZpU?r?x1($) — 2pUKyz2(P), (3.2) 


where vo is a constant. The functions /(¢), 2(d), 71(¢) and y2(¢) are 
found to satisfy the following set of ordinary differential equations: 


POP a ea oe, (3.3) 
%1= 0, (3.4) 
L2= be Ieazte 2) ; (3.5) 
U} h” , 
——h” + 2fh’ = 0, (3.6) 


where S = B2/4supU2. 
The boundary conditions are taken to be 
/(0) = F(0) = 9, f'(4) > 1 as $ +00, (3.7) 
h(0) = 7, h(d) >y asda, (3.8) 
where y and 7 are constants. Equation (3.4) shows that y; is a 
constant, and we see from (3.3) that this must be 
m1 = 1+ SpA; (3.9) 
thus (3.3) becomes 
14 
reais + 2ff! =f? 4-1 = S(h?) = y?). (3.10) 
A suitable solution for the solid region is 


B’g = B'r(1 +..D’2), (3.11) 
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representing an axial current; in order to satisfy the condition 
that the tangential component of the magnetic field should be 
continuous, tr = B’u/Bu’. Note that if y = 7, then A(¢) is constant 
everywhere and (3.10) reduces to the Homann equation on putting 
«x =v. The magnetic field has no effect on the motion in this case, 
whatever the value of S, since on convection a tube of force replaces 
and is replaced by an exactly equivalent tube of force. Convective 
interchange of tubes of force in this manner is important in certain 
geophysical problems 1’), and it has also been discussed by 
Michael }8) in relation to a stability problem. 

As in the previous section, we use F, H, & to denote /, h, 6 when x 
is chosen to be vy, and F, Y, # when « is chosen to be 7. 

Approximate solutions valid for large diffusion 
numbers. For large values of the diffusion number the fluid may 
be considered to be perfectly conducting everywhere outside a thin 
current layer situated next to the boundary. The solution for the 
flow outside this current layer is found by putting « =» and 
neglecting the term H”/e in (3.6) which becomes 


H’(é) =0. (3.12) 


Hence the magnetic field is given by H(é) = constant = y, and 
(3.10) reduces to the Homann equation for F(&). 
In the current layer, it is sufficiently accurate to take F(#) ~ pd? 
where 8 = 0.656/e?; equation (3.6) therefore reduces to 
He?  Q2BPA =A} (Salis) 
which has the solution 
o co 
H (8) = 2 + (y — 7) (f e778 dx) |( fe "8 dx). (3.14) 
0 0 
As was found in the previous section, the ratio of the thickness of 
the current layer to that of the viscous boundary layer is O(e~*) 
when e> 1. 
Approximate solutions valid for small diffusion 
numbers. Proceeding as before, we consider first the current layer 
which envelopes the viscous boundary layer when ¢ <1. Taking 
x = y and neglecting the term eF” in (3.10), we obtain the following 
equations for the current layer: 


LDF. Ure F'2 — t= to — v2), (3.15) 
KH! + IFH' =0. (3.16) 
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In this case the solution in ascending powers of S yields 
F (9) = 09 and Ho(9) = 7 + (y — 7) erf #. Numerical integration 
of these equations for finite values of S can be carried out by the 
method described in the last section, but this is not necessary if 
one is interested only in the viscous boundary layer. 

For the solution in the viscous boundary layer, we put « = » 
and neglect the term ef H’ in equation (3.6). Thus H”(&) = 0 and 
so H(é) =7 to a good approximation in this region. Substituting 
this result in (3.10) we obtain 


BO ee (3.17) 


which reduces to the Homann equation (2.37) on making the 
transformation 


Ele) — LOG) Aaa & (3.18) 
where 
L4 = 1 + S(y? — 7). (3.19) 


From (3.17) we see that F'’(&) + L2 as € + oo, and so the boundary 
conditions on @(A) remain 0(0) = 0’(0) = 0, O’(co) = 1. It appears 
that in cases where 7 > y a sufficiently large value of S causes the 
solution to break down, since L4 becomes negative. In particular, 
if y = 0 and 7 = 1, the breakdown occurs when S > 1, indicating 
that the kinetic energy density of the flow is less than the field 
energy density, and so the fluid cannot convect the current away 
from the stagnation point. It should be noted that there is no 
possibility of reversed flow in this instance as F’(oo) = L?2 becomes 
imaginary at the point of breakdown, and this cannot be altered 
by a simple change of sign in the velocity. 

This result may be of interest to experiments since it shows that 
very strong magnetohydrodynamic boundary layer effects can be 
demonstrated despite the small values of ¢ obtainable in the 
laboratory. The author doubts that a flow could be made to ‘choke 
up’ in practice as an instability would probably occur causing axial 
symmetry to be lost; nevertheless, the trend towards choking could 
be demonstrated for values of S slightly less than unity. 
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THE FLOW OF A CONDUGTING -BLUID PAST: 
MAGNETIZED CYLINDER AT HIGH MAGNETIC 
REYNOLDS NUMBERS 


by L. M. HOCKING 


University College, London, England 


Summary 


A conducting, inviscid fluid flowing past a magnetized cylinder changes 
the magnetic field outside and inside the cylinder. This field is found for a 
cylinder magnetized as a dipole, with the restrictions that (a) the magnetic 
field is weak and (b) the magnetic Reynolds number is large. The magnetic 
field in turn alters the velocity field by the generation of vorticity, and the 
new velocity field and the resulting drag are found when the direction of 
the dipole is parallel to the free stream. 


§ 1. Introduction and formulation. The flow of an inviscid, 
incompressible, conducting fluid past a magnetized cylinder is 
considered in this paper for cases where the ratio of the magnetic 
to the dynamic pressure is small, but the magnetic Reynolds number 
is large. This implies that the magnetic field is considerably dis- 
torted from its value when the fluid is non-conducting, but the 
velocity field is only slightly altered. The magnetic Reynolds 
number is defined by Ry = LUwo, where L and U are a typical 
length and velocity respectively and w is the permeability and o 
the conductivity of the fluid, so that large values of Ry can occur, 
for example, in astrophysical problems, when the large values of L 
will outweigh the small values of o. The natural problem to con- 
sider for astrophysical application would be the flow past a sphere, 
but this is a considerably more difficult problem than the flow past 
a cylinder, and the methods used here cannot be adapted to the 
spherical problem, even when the magnetic field is aligned to 
maintain axial symmetry. 


When Ry is large, the magnetic field tends to be aligned with the 
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velocity, so that a magnetic line emerging from the cylinder is 
swept round the cylinder by the fluid, and the magnetic field is 
confined to a boundary layer and a narrow wake stretching down- 
stream from the rear stagnation point. This magnetic field exerts 
a non-conservative body-force on the fluid, producing vorticity 
and hence modifying the velocity field. It should be made clear that 
the basic velocity field is the simple potential flow with no separation 
of the flow from the cylinder. Obviously, a model which included 
separation would be more realistic but also much more complicated. 
Similar problems to that considered in this paper, but for small 
values of Ry, have been solved by Ludford and Murray 4) and 
Murray and Chi ”) fora sphere and a cylinder, respectively. In the 
present paper, a method is devised for finding the magnetic field 
for cylinders of arbitrary cross-section and for any value of Ry, 
and the particular case of large Ry anda circular cylinder magnetized 
as a dipole is considered. When the dipole is aligned with the free 
stream, the new velocity field and the drag are found. 
The equations satisfied by the velocity and magnetic fields are, 

in non-dimensional form, 

curlg X q = —V(p + 39?) + Pcurl H x A, 

curl H = Kyg * H, 

div H see iy  ==_0, 
where f = wh?/pU2, L, U and h are the units of length, velocity 
and magnetic field, and p is the density of the fluid. As in Murray 
and Chi 2) the electric field has been made zero, which is consistent 
if no current flows at infinity. The parameter f is the ratio of the 
magnetic and dynamic pressure and is supposed to be small so that 
an expansion in powers of 6 can be made. With q = qo + fqi + ---, 
b=f0+ Poi... and H= H+ 6M-+ ..., where qo and fo 
are the velocity and pressure in the absence of any magnetic field, 
the equations for Ho and qi are 


curl Hyp = Rudo X Mo, (1) 

div Hy = 0, (2) 

curl qi X qo = —VPi + furl Ho x Mo, (3) 
div:gi = 0, (4) 


| where P; = 1+ Qo°q1.- 
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It is convenient to introduce a potential A defined by 


0A WH 0A 

H, aa yobs at hw os , 
where and @are polar coordinates. It was pointed out by Michael?) 
that the equation (1) in terms of A, which he derived for a similar 
problem with an externally applied magnetic field, is identical 
with that for temperature in the convection of heat from a cylinder 
by an inviscid stream of fluid. It is well-known that this equation 
can be put in a simple form by using the velocity potential and 
stream function of the irrotational flow as coordinates. For a 
circular cylinder, if 


gy = cos 6 (v2 + 1)/27, y = sin 8 (7? — 1)/2r7, 


so that 2m and 2y are the velocity potential and stream function, 
the equation becomes 

02A 02A OR oA 

pt ape OM a F 
and the boundary of the cylinder in the (9, w)-plane is |p| < 1, 
y = 0. A similar transformation will yield the same equation and 
the same boundary for any cylinder. 

A method of solving this equation in the heat convection problem 
was devised by Piercy and Winny 4) and is used in this paper. The 
method consists of taking a distribution of fundamental singular 
solutions of (5) on the cylinder and using the boundary conditions 
to determine the density of the distribution. The fundamental 
solutions required are 


exp(Raup) Kol Rup? + y?)4), fe (expla) Kol[Ru(y? + y?)*J}, 


which will be called a ‘pole’ and a ‘dipole’ respectively, since they 
reduce to those forms when Ry = 0. The direction of the ‘dipole’ 
is parallel to the p-axis. If a distribution of ‘poles’ only is taken, the 
value of the magnetic potential at any point is 


1 : 
A =s u(é) exp[Ru(p — &)] KofRul(p — £)? + yA} dé, (6) 


where u(é) is the density of the distribution at the point g = é, 
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y = 0 on the cylinder. Similarly, a distribution of ‘dipoles’ only 
gives 
1 


) 
A = | o(8) expikul(p — 8) Fy KotRanlle — 8)? + wR} ae, (2) 
= 

with v(&) as the density of the distribution. 

The ‘pole’ distribution (6) makes A an even function of y. The 
value of 0A/éy on p = 0 is O if |g| > 1, but if |p| < 1, 

if 

cA) = tim fu(é) expiR : 
oa = hm: | “(é) ex ; é eee 
OY |y-0 v0 BANE i oe) eat wall 


=! 


» K’otRaal(p — €)? + w?]}} dé. 
Expanding the integrand, the only singular part comes from the 
neighbourhood of the point & = q, so that 


1 


= lim — yu(g) 


y=0 y—0 


aA 
op 


dé ie | —zu(p) 1f yO, 
g—é)j2+y? | +au(9) if y+O0-. 


(8) 


Similarly, for the ‘dipole’ distribution (7) 


A|,9 = F (9), (9) 
and 
oA o2 
i. | 8 exrale — 8) yp HolRadle— a+ NA 
= (8)(2Rue = © )lexpiRarl—8)1K ol Rarip E108, (10) 
be p ay? 
using (5). 


As in Murray and Chi’s paper the magnetic field is supposed 
to originate inside the cylinder as a dipole field. If the direction of 
the dipole is perpendicular to the direction of the free stream, A 
must be an even function of y by symmetry, and the solution can 
be found by a distribution of ‘poles’ only. When the dipole is 
parallel to the free stream, A is odd in y, and the ‘dipole’ distribution 
(7) must be used. Since (5) is linear in A, a dipole at any angle 


234 L. M. HOCKING 


requires a linear combination of the two special cases, so that only 
these need be considered. 

The conductivity of the cylinder is finite, so that the field inside 
the cylinder is modified by the motion of the conducting fluid 
outside. If A’ determines the field inside, A’ is a harmonic function 
which contains a singular term, giving the dipole field, as well as 
terms finite inside the cylinder. The boundary conditions which 
link A and A’ are the continuity of tangential magnetic field and 
normal magnetic induction. Thus, in the (7, 0) plane, on 7 = 1 

0A 0A’ eS 

=p eral pAs—w Ad, (11) 
where ww and yw’ are the permeabilities of the fluid and the cylinder 
respectively. Also, on, 7 = 15 


0A 0A op 0A dp J 0A 


§ 2. Dipole pointing upstream. When the dipole in the cylinder 
is pointing upstream, the appropriate form for the field inside is 


sin 0 


B + ¥ @nr™ sin nO, (12) 
n=1 


A 


and for the field outside, A is given by (7). The boundary conditions 
on the upper surface of the cylinder give 
; , 
é 02 
Joe (28 Z — ) exptRnly — £)) KolRalp — 8) a = 
i es ms 


sinO a |,’ Us) 
hoes 


with 0 <6 <a and gy =cos§. The boundary conditions on the 
lower surface are also satisfied since A and A’ are even functions of y. 
If Ry is small, the solution can be found by expanding in powers 
of Ry. However, this is the case considered by Murray and Chi, 
and the solution obtained by the present method agrees with their 
solution so the details will not be given here. It would not be diffi- 
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cult to find further terms in the Ry expansion, but the main 
features of the magnetic field for small Ry are determined by the 
term already found, so that it is more interesting to consider large 
values of Ry. 

The asymptotic value of the Bessel function for large Ry and 
lp — &| not small is 


Ko(Rarip—&) = (52) (ip Ei) texp(—Rarip—E {1 + (Ra) 


If the contribution from the small neighbourhood of the point = 9, 
in which Ry |p — &| is not large, is neglected, the integrand in(13) is 


(21Ra)? 


(89 — II 
ip 


if € < m and is exponentially small if § > gy, so that (13) becomes 


2nR ao dé = wate (15) 
(2x7Rm)* 0b es, arbi ay = 
To solve this ree and (14), write 
aust sin 
1 
v= Rub U1 “+ Bore 
so that the leading terms are given by 
"nlé) 
re) V1 
22)* | de 2, (16) 
op Se —® 
1 
TU, 
Malate — yp vue). (12) 


The integral equation for v; can be written 


[28 a= (2) +n. 
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where k is some constant, and the solution of this equation is 


no 2(2) See 2 (ery, 0 


where & has been chosen so that v1 is finite at & = —1. The value 
of A’; on the cylinder is given by (17) as 


Ata er 


4 
P cos 40,0< 0 =< a, (19) 
mo 


and since A’ is an odd harmonic function, the value of the potential 
inside the cylinder is easily found to be 
sin 6 32u = 


A’ = — vy sin 9 — —--- 2S 
Y | mr’ (a7 )* a 4n2 — 


; vsin nO. (20) 


The value of A so far obtained has a discontinuity at the rear 
stagnation point, gy = 1. Mathematically, this has appeared because 
the contribution of the neighbourhood of the point = @ to the 
integral (13) cannot be neglected when gy = 1. The magnetic field 
is shown by (20) to be almost confined to the interior of the cylinder, 
but there are magnetic lines, for small values of A of order Ry“, 
which leave the cylinder at all points, re-entering the cylinder at 
the rear stagnation point. If the neglected part of (13) were included, 
it would show that the re-entry took place over a region near the 
stagnation point with a width of order Ry“. 

The field outside the cylinder is given by (7) and (18). Using the 
asymptotic expansion of the Bessel function to the same order as 
used in satisfying the boundary conditions, A = 0 if pg < —1, “.e. 
upstream from the cylinder and, if |p| < 1, 


alg! [(p — —F 2 4 vo exp[ — Ruy?/2(p — &)] dé. (21) 


The magnetic field is negligible outside the boundary layer defined 
by 
Ruy? = 2a(p + 1) (22) 


where A is some constant such that exp / is large compared with 1. 
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In the effective part of the integral in (21), Ru(p — €) is large 
and »/(y — &) is small. Defining new variables by 


SA: Ruy? 
Ca, gf (23) 
2p’ 2(p + 1) 


and neglecting y in comparison with (g — &) in the denominator, 


af ate +) {( ) 
A =—_ — | =f $e: =o Sa 
2 ae l iD exp(— ¢*) d¢. (24) 


£1 


After some manipulation, this can be written 


co 


Dlw27 3 
A= — =| ED | exp(— 212) — 24 fexp(—e2) atl. (25) 


a1 


I 


On the cylinder z; = O, and (25) gives the same value of A as can 
be found from (19). 
Downstream from the cylinder, 


(1 + €)ty 


JU o — &)2 + ye} exp[— Ruy?/2(p — &)] dé (26) 


with m > 1, so that A is negligible outside a wake defined by the 
same expression as the boundary layer (22). In the wake, by the 
same methods as those used above, the value of A is 


7 qe en) Ti a) exp(— ¢2) dé, (27) 


where 2; is defined in (23) and 


Ruy? 
2(p — 1) © 


At points in the wake sufficiently far from the cylinder for g — & 
to be replaced by its average value in (26) 


a 8 
A= Ses vt exp(— Rarv?/2e), 


222 => 
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‘ A A : 1 Senet 
or, since in this region gy ~ 47 cos 6 = 4x and y= 9 sin 0 = 9y, 


16y 
A = — —— exp(— Ruy?/4x). (28) 
3x? 


Traversing the positive side of the wake, keeping x fixed, A de- 
creases from 0 on the axis to a minimum value 


16 aa. : 
male sees 
and then tends rapidly to 0. As Ry increases, each magnetic line 
moves towards the cylinder and narrows until A = 0 everywhere 
outside the cylinder. The limiting value of the magnetic field as 
Ry tends to infinity is zero except on the cylinder, where there is 
a tangential component of magnitude 2sin @ directed towards the 
rear stagnation point, and on the downstream axis, where there is 
a field of magnitude 16/3? for large x, directed towards the cylinder. 

A sketch of the magnetic lines for large Ry is shown in fig. 1. 


Fig. 1. The magnetic lines when the dipole is aligned with the stream. The 
full lines are calculated from equations (25) and (28) for Rag = 100, and the 
broken lines are sketched to connect the calculated values. 


§ 3. Dipole perpendicular to the free stream. When the dipole is 
at right angles to the direction of the free stream, the appropriate 
form for the field inside is 


ae cos 6 


+ ¥ anr®™ cos nO, 
n=1 
and the field outside is given by (6). The boundary conditions 
on the upper surface of the cylinder give 
| me Be 


— 1Uu(y) = — — 
7) sin@ Or \,04 


P (29) 


[u&) exp(Ru(p — 6] KolRu lp — 1) ae = wa’ 


=i 


(30) 


r=1° 
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As before, the solution for small Ry, found by an expansion in 
powers of Ry, yields the same results as were found by Murray and 
Chi. For large values of Ry the asymptotic expansion of the 
Bessel function replaces (30) by 


- 
( Mi 4 )| Aas ; dé ay, nt (31) 
2Ku ip fy le 
=! 
Writing 
cos 
ye le a apie 
1 
u(p) = u(y) + Ryr vy) + 
M” 
the leading term in (29) gives 
2 
m(9) = a, (32) 
(1 — ?)? 


and the value of A on the cylinder is 


o se c= ae Be i 


This integral can be expressed in terms of K and E, the first and 

second complete elliptic integrals, giving 
2 

A = — ——— [K(cos 36) — 2E(cos $6)], (34) 

where cos § = g. Thus on the cylinder the potential varies from 

(z/Rm)! at the front stagnation point, 9 =z, to O at the point 

§ ~ 49° and is negative for smaller values of 0. The expression 

(34) tends to infinity at the rear stagnation point. The correct value 
at this point, 6 = 0, is given by (30) as 


hayes 2 Fat vy exptRa(! — 6)] KolRan(t — 2)) 48, (85) 


and: the divergent integral arises when the asymptotic value of 
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the Bessel function is used at € = 1. To evaluate (35) more accurate- 
ly, the integral can be divided into two parts, in one of which the 
asymptotic expansion can be used. If it is supposed that the value 
of Ko(x) is given by the first term of the asymptotic expansion to 
sufficient accuracy when x > é: 
1—e/Ru 1 
4=—-(35) Jaceroce alo 
~  \aku/ J (1-0 — 2b od (1 — 8) 
=i 1—e/Ru 


- exp [Ru(1 — )] Ko[Rm(1 — §)] dé = 


l 
— RENE oe Ru — O(Ru), 

on evaluating the first integral, the second being of order Rag7?}. 
The value of the potential on the cylinder thus tends to 0 as Ryy + co 
at all points, but it becomes relatively more concentrated near the 
rear stagnation point. 

As before, A is exponentially small except within the boundary 
layer and the wake, defined by (22). In the boundary layer, 


DNs é 
€ (=) | Be er Ceili aie 
-1 

and the value in the wake is given by the same expression with the 
upper limit of integration 1. The forms obtained for these integrals 
by making the transformations (23) are not particularly simple 
and so are omitted. The value of the potential in the wake at a large 
distance downstream is found to be 


WA + 

Ass — feat x—* exp(— Ryuy2/4x) (37) 
Ru 

in the same way as (28) was found. 

In the limiting case, as Ry tends to infinity, there is a tangential 
component of magnetic field, Hg = —2 cos 0, on the cylinder. In 
the wake, the y-component of the magnetic field tends to 0, but 
the x-component ' 

0A 


Ay =r 
--—=( 


aR $ 
7) yx exp(— Ruy?/4z) 
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has a maximum value (z/2e)! «2 at y = (2x/Ry)* and a correspond- 
ing minimum at y = — (2x/Ry)!. In the limit, the magnetic field 
in the wake lies along the axis and is directed away from and towards 
the cylinder on opposite sides of the axis. A sketch of the magnetic 
lines for large Ry is given in fig. 2. 


Fig. 2. The magnetic lines when the dipole is perpendicular to the stream. 
(N.B. This is a rough sketch intended to give qualitative information only). 


§ 4. The modified velocity field and the drag. The magnetic field 
in the fluid produces a non-conservative body force which generates 
vorticity and so modifies the velocity field. The equation determin- 
ing this modification is (3), which leads to 


pee Ryne. 38 

w 2 bp + KuMligAy (38) 
ligr 

ee ENR Ps (39) 
2 Op 


where w is the vorticity, perpendicular to the plane of motion, 
pi and q; are the additional pressure and velocity terms, P} = ~1 + 
+ qo'qi, and Hg and Hy are the g- and y-components of the 
magnetic field. Since » must vanish at a large distance upstream, 
i.e. at p = — oo, the terms in P; can be eliminated giving 


P 
0 
id ~ [ageap, 40 
wef yp AP (40) 


siice 40 when go =< — I. 

The vorticity is negligible except in the foutialeesey: layer and the 
wake, a narrow region on either side of the line y = 0. This surface 
distribution of vorticity can be replaced by a line distribution on 
_y =0, the strength at any point being determined by integrating 
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w across the region, 7.e. from y = 0 to y = co on the upper surface 
of the cylinder, from y = 0 to y= — oo on the lower surface, 
and from y = — oo to y = co in the wake. From this distribution 
of vortex filaments, the velocity can be found and an irrotational 
velocity field can then be added to satisfy the boundary conditions 
on the cylinder and at infinity. 

For the dipole pointing upstream, the magnetic field is determined 
by (24). The vorticity density on the upper surface is given by 


@ =fo dy = Ruf Hy dy + (Ror [Hy deli=0 


using (40), and substituting the values found from (24), 
~ 2 l 2 / 
Den Verupirermtas opore a) (41) 


On the lower surface, the sign is reversed, and since the vorticity 
has opposite signs on either side of the downstream axis, the 
vorticity density along the axis is zero. The velocity field produced 
by this vorticity distribution can be found in the usual way by 
summing the velocities produced at any point by each vortex 
filament. On the cylinder, the velocity produced at the point 6 
by the filaments at 6 == + 6; is directed radially outwards and is 
of magnitude 


&(cos 61) 


5. [cot $(01 — 6) + cot 401 + )], 
JU 


so that the velocity field produced by the vorticity has only a 
radial component on the cylinder, 
me Es | (cos 61) sin 61 


phe 27 cos'@; "cos 0 
0 


dy. (42) 


Since the total strength of the vorticity is zero, the velocity is 
O(1/r?) at large distances from the cylinder. The irrotational field 
to be added to this must vanish at infinity and make the total radial 
component zero on the eyvidnder: It can be derived from a potential 
of the form 


© ay, cos nO 
KT ee ae ee ot Re 


n=1 aS 
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and the boundary condition dy/ér = — 9, on r= 1 gives 
co 1 [ (cos 63) sin 01 
X Nan Cos nO = | d6;, 43 
4x4 i 27 cos 6; — cos 0 ; Ge} 


0 


which determines the coefficients @p. 


The pressure on the cylinder is found by integrating (39), which 
gives 


é, ] 
di + 2sin 8 } nay sin nd = — Rar | Hy? dp = — — (3+ 2p — 9). 
4 
=I 


The drag on the cylinder is produced by the pressure and the 
Maxwell stress. The pressure drag can be written pU2LDp, where 
Dp is a non-dimensional drag coefficient 


Dp = — 8 me cos 6 dé = B(2 + 2na»). (44) 
0 


The coefficient az is found from (43) to be 


@(cos 61) sin 6; cos 26 4 
ao —— dé do; a=) WES 
222 cos 6; — cos 6 322 


so that 

8 
pees (45) 

om 

The Maxwell stress produces a drag with coefficient 

Qa 
Du = B / [4(A,? — Ho?) cos 0 — H,Hg sin 6] dé, 
0 


and substituting the value of the magnetic field on the cylinder, 


pees (A) e ie (46) 
15 ww’ Ru 

The pressure drag is of larger order than the Maxwell drag and is 

independent of Ry, so that (45) is the value of the drag when Ry 

is infinite. This drag is produced by the magnetic forces, not directly 

by means of the Maxwell stress, but through the production of 

vorticity and the consequent modification of the pressure. It is 
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interesting to note that the drag on the cylinder for infinite Fy, is 
proportional to BpU?2L = wh?L and so is independent of the stream 
velocity. This does not mean that the force is non-zero when the 
fluid is at rest, since the analysis is based on the supposition that 
uh?/pU2 is small. 

When the dipole is perpendicular to the free stream, the more 
complicated form for the magnetic field and the difficulty of re- 
moving the infinity at the rear stagnation point make the determina- 
tion of the velocity field and the pressure too difficult to be ac- 
complished by analytic methods. The Maxwell drag can be found 
using (29) and (34), giving a drag coefficient 


Dy = 4.158 (4 = ) Rar i, 
lu 


but the pressure drag may be expected to be of higher order than 
this, as in the previous case. 

I would like to express my indebtedness to Dr. J. D. Murray 
for his help and encouragement during the preparation of this paper. 
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One of the 5 types of X-ray microscopes is the projection 


microscope, originally proposed by Von Ardenne in 1939 
and realized by Cosslett and Nixon in 1951. 
With the aid of an ultra-fine-focus X-ray tube with a spot 


diameter of 1-0.1 micron, an enlarged X-ray shadow 


image is projected on a screen or film. Due to negligible 


scattering and refraction clear absorption images can be 


obtained. A special feature is the great depth of focus. 


Thus ideal stereo images can be made, showing the three 
dimensional mass absorption distribution of the speci- 


men. This will be of great importance for both research 


(medical, biological, metallurgical, criminological, etc.) 
and technology (rubber-, textile-, paper-industry etc.). 
In “Microprojecton with X-rays” the author describes 
his experiences in constructing and applying the pro- 
jection microscope in practice, covering five years of 
research in this field. His book will serve as a guide to 
those who are intending to explore the possibilities of this 
type of microscope. 
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